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(57) Abstract 

Transmission of a sequence 
of data samples includes transmis- 
sion of a sandwich preamble ("san- 
damble") that has a first pream- 
ble that precedes the sequence of 
data samples, and a second pre- 
amble that follows the sequence 
of data samples. The sandamble 
may include several guard regions, 
including an extended guard re- 
gion. Various repetitions of data 
in the sandamble make it useful 
for aiding time and frequency syn- 
chronization of the received sig- 
nal when it has been subjected 
to a dispersive channel. In an- 
other aspect, a coarse estimate of 
a frame start position of a re- 
ceived sequence of desired data 
samples in a received signal se- 
quence may be improved by using 
the coarse estimate to generate fre- 
quency-domain received samples. 
A frequency-domain correlation is 
then determined between the fre- 
quency-domain received samples 
and noiseless samples. The im- 
proved estimate of the starting position of the sequence of desired samples in the received signal sequence may then be generated based on 
an argument of the frequency-domain correlation between the frequency-domain received samples and the noiseless samples. The noiseless 
samples may be generated by demodulating symbols contained in the received signal sequence. In alternative embodiments, the noiseless 
samples may be generated from pilot symbols that have been multiplexed into the sequence of received samples. 
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BURST CARRIER FREQUENCY SYNCHRONIZATION AND 
ITERATIVE FREQUENCY-DOMAIN FRAME SYNCHRONIZATION 

FOR OFDM 

CROSS-REFERENCE TO COPENDING APPLICATION 

This application claims the benefit of U.S. Provisional Application 
No. 60/088,438, filed June 8, 1998, the entire disclosure of which is hereby 
incorporated herein by reference. 

BACKGROUND 

The present invention generally relates to radio communications 
systems, and more particularly to frame and frequency synchronization of bursts 
received via a dispersive channel. 

Performing frame and frequency synchronization of transmissions 
that have been received over unknown frequency-selective channels (i.e., 
dispersive channels that cause Inter-Symbol Interference (ISI)) is a problem that 
calls for different solutions than are used when the transmission takes place over a 
non-selective channel. For time synchronization on frequency-flat (i.e., non- 
selective) channels, time synchronization via peak detection in the receiver is 
usually performed by a correlation filter that detects a specific correlation sequence 
that is inserted by the transmitter. But this procedure will not produce a distinct 
peak for transmission over unknown frequency-selective (ISI) channels. In other 
words, the formerly good (optimized) correlation properties are destroyed by the 
convolution of the originally transmitted signal (including the specific correlation 
sequence) with the unknown channel impulse response. Apart from this problem, 
this preamble is not useful for performing carrier frequency synchronization. 
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For continuous transmission, such as in broadcast applications, the 
receiver can average the synchronization parameters like frame start position and 
frequency offset over several preambles in order to obtain a very accurate and 
reliable frame and carrier frequency synchronization result. 

A severe problem in other systems, such as in wireless 
Asynchronous Transfer Mode (ATM) scenarios, arises as a result of the packet- 
oriented transmission and the mostly non-continuous traffic. This requires mostly 
burst synchronization schemes that allow a reliable single-shot frame- and carrier 
frequency synchronization to be performed. 

For frame and carrier frequency synchronization of spontaneous 
transmissions taking place over unknown frequency-selective fading channels, a 
special preamble structure has been proposed which consists of some channel 
symbol sequence that is repeated one or more times, so that periodicity is 
introduced into the transmitted signal. This is described, for example, in Pierre R. 
Chevillat, Dietrich Maiwald, and Gottfried Ungerboeck, "Rapid Training of a 
Voiceband Data-Modem Receiver Employing an Equalizer with Fractional- r 
Spaced Coefficients", IEEE Transactions on Communications , vol. 35, no. 9, 
pp.869-876, 1987 (henceforth "[CMU87]"); Stefan A. Fechtel and Heinrich Meyr, 
"Fast Frame Synchronization, Frequency Offset Estimation and Channel 
Acquisition for Spontaneous Transmission over Unknown Frequency-Selective 
Radio Channels", Proceedings of the International Symposium on Personal. Indoor 
and Mobile Radio Communications (PIMRC^ . PP. 229-233, Yokohama, Japan, 
1993 (henceforth "[FM93]"); Stefan A. Fechtel and Heinrich Meyr, "Improved 
Frame Synchronization for Spontaneous Packet Transmission over Frequency- 
Selective Radio Channels", Proceedings of the International Symposium on 
Personal, Indoor and Mobile Radio Communications (PIMRC'94) T pages 353-357, 
The Hague, Netherlands, 1994 (henceforth "[FM94]"); and Uwe Lambrette, 
Michael Speth, and Heinrich Meyr, "OFDM Burst Frequency Synchronization by 
Single Carrier Training Data", IEEE Communications Letters , vol. 1, no. 2, pp. 
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46-48, 1997 (henceforth "[LSM97]"). This type of preamble is referred to herein 
as a "repetition preamble" . Examples of these conventional repetition preamble 
structures are given in FIGS, la and lb. In the conventional repetition preamble 
depicted in FIG. la, the transmitted signal is replicated in the regions designated A 
5 and B, and the guard region, G, is usually a copy of the rightmost part of the A 

region. The conventional repetition preamble depicted in FIG. lb is similar, but 
here the transmitted signal is replicated more than once; that is, the signal is 
identical in each of the regions A, B and C, with the guard region, G, again 
usually being a copy of the rightmost part of the A region. In each case, the 

10 replicated regions (i.e., A and B, or A, B and C) are contiguous to one another. 

The data to be transmitted in the frame is arranged so that it follows all of the 
replicated preamble regions. 

After convolution of the periodic signal part with the (finite) 
impulse response of the unknown frequency-selective (ISI) channel, the received 

15 signal in the regions A and B (or A, B and C for the case of FIG. lb) will again 

exhibit some similarity, assuming that the preamble part G is chosen to be 
sufficiently long. This is true, even though the shape of the received signal in 
these regions can be completely different from the transmitted one due to the 
frequency selectivity (time dispersivity) of the channel. The only difference 

20 between the received signal in regions A and B (and B and C) will be a phase shift 

that is proportional to the carrier frequency offset. 

Thus, the receiver can detect the correct starting position of the 
received signal by processing the received samples, including performing a signal 
correlation, given that the preamble samples are spaced apart by the discrete 

25 periodicity interval Icq. This is described in Jan-Jaap van de Beek, Magnus 

Sandell, Mikael Isaksson and Per Ola Borjesson, "Low-Complex Frame 
Synchronization in OFDM Systems", Proceedings of the International Conference 
on Universal Personal Communication (ICUPC'95) . pp. 982-986, Tokyo, Japan, 
1995 (henceforth "[vdBSIB95]"); Magnus Sandell, Jan-Jaap van de Beek, and Per 
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Ola Borjesson, "Timing and Frequency Synchronization in OFDM Systems Using 
the Cyclic Prefix", Proceedings of the International Symposium on 
Synchronization , pp. 16-19, Essen, Germany, 1995 (henceforth "[SvdBB95]"); 
Timothy M. Schmidl and Donald C. Cox, "Low-Overhead, Low-Complexity 
[Burst] Synchronization for OFDM", Proceedings of the international Conference 
on Communications QCC'96^ pp. 1301-1306, Dallas, Texas, USA, 1996 
(henceforth "[SC96]"); and [LSM97]. 

To utilize the conventional repetition preambles of FIGS, la and lb, 
the signal parts in regions A and B (and B and C and eventually A and C) are 
processed to obtain the desired synchronization parameters time and carrier 
frequency. The exploitable periodicity interval(s) is (are) illustrated by the lines 
101, 103 and 105. A minimum (or maximum) value of the timing metric occurs 
not only at the correct time position, but also in a wide range around it. Hence, 
these conventional preambles suffer from the same ambiguity problem as the 
correlation-sequence technique in ISI channels. (In this context, ambiguity should 
be understood as a blurred extreme point of the timing metric.) If additionally a 
high noise power is present at the receiver input, the probability of an error in time 
synchronization is high, resulting in a very high variance of the timing estimate. It 
should be noted that, as explained in [SvdBB95] and [SC96], the argument of the 
correlation result between A and B at the correct timing instant offers an estimate 
for the frequency offset. Thus, it can be exploited for frequency synchronization 
purposes. With respect to this point, the repetition preamble offers at least one 
advantage over the correlation sequence preamble. 

SUMMARY 

It is therefore an object of the present invention to provide 
apparatuses and methods for performing frequency synchronization of received 
signals. 
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It is a further object of the present invention to provide apparatuses 
and methods for performing frame synchronization of received signals. 

The foregoing and other objects are achieved in methods and 
apparatuses for transmitting and receiving a sequence of data samples. In 
accordance with one aspect of the invention, a sequence of data samples is 
transmitted by initially transmitting a first preamble comprising a sequence of 
preamble samples, then transmitting the sequence of data samples, and then 
subsequently transmitting a second preamble comprising the sequence of preamble 
samples, whereby the sequence of data samples is transmitted after the step of 
initially transmitting the first preamble, and before the step of subsequently 
transmitting the second preamble. 

In another aspect of the invention, the first preamble comprises a 
first guard region and a symbol sequence; the first guard region comprises a first 
subset of the symbol sequence; the second preamble comprises a second guard 
region and the symbol sequence; and the second guard region comprises a second 
subset of the symbol sequence. 

In yet another aspect of the invention, the first subset of the symbol 
sequence may be a larger subset of the symbol sequence than the second subset of 
the symbol sequence. In alternative embodiments, the second subset of the symbol 
sequence may be a larger subset of the symbol sequence than the first subset of the 
symbol sequence. 

In still another aspect of the invention, a sequence of desired data 
samples is received by receiving a signal sequence that includes a first preamble 
followed by the sequence of desired data samples followed by a second preamble, 
wherein the first preamble comprises a symbol sequence, and the second preamble 
comprises the symbol sequence. A plurality of locations of the first preamble in 
the received plurality of data samples are then hypothesized. For each of the 
hypothesized locations, a hypothesized first preamble and a corresponding 
hypothesized second preamble are determined. For each of the hypothesized 
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locations, a correlation between the hypothesized first preamble and the 
corresponding hypothesized second preamble is determined. The correlations are 
used to determine locations of a most likely hypothesized first preamble and of a 
most likely second preamble. 

In yet another aspect of the invention, the locations of the most 
likely hypothesized first and second preambles are used to determine a first 
frequency offset of the received signal sequence. 

In still another aspect of the invention, the first preamble may 
comprise a first guard region and a symbol sequence; the first guard region may 
comprise a first subset of the symbol sequence and a second subset of the symbol 
sequence; the second preamble may comprise a second guard region and the 
symbol sequence. Furthermore, receiving the sequence of desired data samples 
farther involves using the first frequency offset to make a first frequency 
correction of the received signal sequence. A frequency corrected first guard 
region is then determined in the first frequency corrected received signal sequence, 
and a frequency corrected first preamble is determined in the first frequency 
corrected received signal sequence. The second subset of the symbol sequence 
from the frequency corrected first guard region is then correlated with the second 
subset of the symbol sequence from the frequency corrected first preamble, 
thereby generating a correlation result. The correlation result is then used to 
detect a second frequency offset. 

In yet another aspect of the invention, using the correlation result to 
detect the second frequency offset comprises using a priori knowledge of 
probabilities of possible frequency offset intervals to detect the second frequency 
offset. 

In an alternative embodiment, in which the first preamble comprises 
a first guard region and a symbol sequence; the first guard region comprises a first 
subset of the symbol sequence and a second subset of the symbol sequence; and the 
second preamble comprises a second guard region and the symbol sequence, 
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receiving the sequence of desired data symbols further includes correlating the 
second subset of the symbol sequence from the first guard region with the second 
subset of the symbol sequence from the first preamble, thereby generating a 
correlation result; using the first frequency offset to correct the correlation result; 
5 and using the corrected correlation result to estimate a second frequency offset. 

In still another aspect of the invention, the first preamble comprises 
a first guard region and a symbol sequence; the first guard region comprises a first 
subset of the symbol sequence and a second subset of the symbol sequence; the 
second preamble comprises a second guard region and the symbol sequence. 

1 0 Receiving the sequence of desired data symbols further includes using the first 

frequency offset to make a first frequency correction of the received signal 
sequence; determining a frequency corrected first guard region in the first 
frequency corrected received signal sequence; determining a frequency corrected 
second preamble in the first frequency corrected received signal sequence; 

1 5 correlating the second subset of the symbol sequence from the frequency corrected 

first guard region with the second subset of the symbol sequence from the 
frequency corrected second preamble, thereby generating a correlation result; and 
using the correlation result to detect a second frequency offset. 

In yet another aspect of the invention, using the correlation result to 

20 detect the second frequency offset comprises using a priori knowledge of 

probabilities of possible frequency offset intervals to detect the second frequency 
offset. 

In an alternative embodiment, in which the first preamble comprises 
a first guard region and a symbol sequence; the first guard region comprises a first 
25 subset of the symbol sequence and a second subset of the symbol sequence; and the 

second preamble comprises a second guard region and the symbol sequence, 
receiving the sequence of desired data symbols further includes correlating the 
second subset of the symbol sequence from the first guard region with the second 
subset of the symbol sequence from the second preamble, thereby generating a 



WO 99/65180 



PCT/EP99/03793 



-8< 

correlation result; using the first frequency offset to correct the correlation result; 
and using the corrected correlation result to estimate a second frequency offset. 

In still another aspect of the invention, the first preamble comprises 
a first guard region and a symbol sequence; the first guard region comprises a first 
subset of the symbol sequence and a second subset of the symbol sequence; and the 
second preamble comprises a second guard region and the symbol sequence. Here, 
the technique for receiving the sequence of desired data symbols further includes 
using the first frequency offset to make a first frequency correction of the received 
signal sequence. A frequency corrected first guard region is determined in the 
first frequency corrected received signal sequence, and a frequency corrected first 
preamble is determined in the first frequency corrected received signal sequence. 
Also, a corrected second preamble is determined in the first frequency corrected 
received signal sequence. The second subset of the symbol sequence from the 
frequency corrected first guard region is correlated with the second subset of the 
symbol sequence from the frequency corrected first preamble, thereby generating a 
first correlation result. Also, the second subset of the symbol sequence from the 
frequency corrected first guard region is correlated with the second subset of the 
symbol sequence from the frequency corrected second preamble, thereby 
generating a second correlation result. The first and second correlation results are 
then combined, thereby generating a combined correlation result. The combined 
correlation result is used to detect a second frequency offset. 

In yet another aspect of the invention, the step of using the 
combined correlation result to detect the second frequency offset comprises using a 
priori knowledge of probabilities of possible frequency offset intervals to detect the 
second frequency offset. 

In an alternative embodiment, in which the first preamble comprises 
a first guard region and a symbol sequence; the first guard region comprises a first 
subset of the symbol sequence and a second subset of the symbol sequence; and the 
second preamble comprises a second guard region and the symbol sequence, the 
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technique for receiving the sequence of desired data symbols further includes 
correlating the second subset of the symbol sequence from the first guard region 
with the second subset of the symbol sequence from the first preamble, thereby 
generating a first correlation result; correlating the second subset of the symbol 
5 sequence from the first guard region with the second subset of the symbol 

sequence from the second preamble, thereby generating a second correlation 
result; using the first frequency offset to correct the first correlation result; using 
the first frequency offset to correct the second correlation result; combining the 
first and second corrected correlation results, thereby generating a combined 

10 corrected correlation result; and using the combined corrected correlation result to 

estimate a second frequency offset. 

In still another aspect of the invention, the received signal sequence 
includes the first preamble, followed by a cyclic prefix guard interval, followed by 
the sequence of desired data samples, followed by a cyclic postfix guard interval, 

15 followed by a second preamble. Furthermore, receiving the sequence of desired 

data symbols further comprises using the location of the most likely hypothesized 
first preamble to determine a coarse estimate of a starting position of the sequence 
of desired samples in the received signal sequence. An improved estimate of the 
starting position of the sequence of desired samples in the received signal sequence 

20 is then generated by: using a time domain to frequency domain transformation 

technique and the coarse estimate of the starting position of the sequence of desired 
samples in the received signal sequence to generate frequency-domain received 
samples; determining a frequency-domain correlation between the frequency- 
domain received samples and noiseless samples; and generating the improved 

25 estimate of the starting position of the sequence of desired samples in the received 

signal sequence based on an argument of the frequency-domain correlation 
between the frequency-domain received samples and the noiseless samples. 

In yet another aspect of the invention, the cyclic prefix guard 
interval comprises an ending subset of the sequence of desired data samples; and 
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the cyclic postfix guard interval comprises a beginning subset of the sequence of 
desired data samples. 

In still another aspect of the invention, the number of samples in the 
cyclic prefix guard interval and the number of samples in the cyclic postfix guard 
5 interval are optimally adapted to an expected channel power delay profile. 

In yet another aspect of the invention, receiving the sequence of 
desired data symbols further includes generating a further improved estimate of the 
starting position of the sequence of desired samples in the received signal sequence 
by: using the improved estimate of the starting position of the sequence of desired 

10 samples in the received signal sequence to generate improved frequency-domain 

received samples; determining a second frequency-domain correlation between the 
improved frequency -domain received samples and the noiseless samples; and 
generating the further improved estimate of the starting position of the sequence of 
desired samples in the received signal sequence based on an argument of the 

15 second frequency-domain correlation between the frequency-domain received 

samples and the noiseless samples. 

In still another aspect of the invention, the noiseless samples may be 
generated by demodulating symbols contained in the received signal sequence. 
This may include generating the noiseless samples by demodulating symbols 

20 contained in the sequence of desired samples. In alternative embodiments, the 

noiseless samples may be generated from pilot symbols that have been multiplexed 
into the sequence of received samples. 

In yet another aspect of the invention, receiving a sequence of 
desired data samples is accomplished by receiving a signal sequence that includes a 

25 cyclic prefix guard interval followed by the sequence of desired data samples, 

followed by a cyclic postfix guard interval, a coarse estimate of a starting position 
of the sequence of desired samples in the received signal sequence is determined, 
and an improved estimate of the starting position of the sequence of desired 
samples in the received signal sequence is generated by: using the coarse estimate 
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of the starting position of the sequence of desired samples in the received signal 
sequence to generate frequency-domain received samples; determining a 
frequency-domain correlation between the frequency-domain received samples and 
noiseless samples; and generating the improved estimate of the starting position of 
5 the sequence of desired samples in the received signal sequence based on an 

argument of the frequency-domain correlation between the frequency-domain 
received samples and the noiseless samples. 

In still another aspect of the invention, the cyclic prefix guard 
interval comprises an ending subset of the sequence of desired data samples; and 

10 the cyclic postfix guard interval comprises a beginning subset of the sequence of 

desired data samples. 

In yet another aspect of the invention, the number of samples in the 
cyclic prefix guard interval and the number of samples in the cyclic postfix guard 
interval are optimally adapted to an expected channel power delay profile. 

15 In yet another aspect of the invention, receiving the sequence of 

desired data symbols further includes generating a further improved estimate of the 
starting position of the sequence of desired samples in the received signal sequence 
by: using the improved estimate of the starting position of the sequence of desired 
samples in the received signal sequence to generate improved frequency-domain 

20 received samples; determining a second frequency -domain correlation between the 

improved frequency-domain received samples and the noiseless samples; and 
generating the further improved estimate of the starting position of the sequence of 
desired samples in the received signal sequence based on an argument of the 
second frequency -domain correlation between the frequency-domain received 

25 samples and the noiseless samples. 

In still another aspect of the invention, receiving the sequence of 
desired data symbols further includes using the further improved estimate of the 
starting position of the sequence of desired samples in the received signal sequence 
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to adjust a measurement of frequency offset associated with the sequence of 
desired samples in the received signal sequence. 

In yet another aspect of the invention, the noiseless samples may be 
generated by demodulating symbols contained in the received signal sequence. 
5 This may include generating the noiseless samples by demodulating symbols 

contained in the sequence of desired samples. 

In alternative embodiments, the noiseless samples may be generated 
from pilot symbols that have been multiplexed into the sequence of received 
samples. 

10 In still another aspect of the invention, receiving the sequence of 

desired data symbols further includes using the improved estimate of the starting 
position of the sequence of desired samples in the received signal sequence to 
adjust a measurement of frequency offset associated with the sequence of desired 
samples in the received signal sequence. 

15 In yet another aspect of the invention, receiving the sequence of 

desired data symbols further includes iteratively generating a final improved 
estimate of the starting position of the sequence of desired samples in the received 
signal sequence by performing a sequence of iterative steps for a number of times. 
Each sequence of iterative steps comprises using a previously generated estimate of 

20 the starting position of the sequence of desired samples in the received signal 

sequence to generate current frequency-domain received samples; determining a 
current frequency-domain correlation between the current frequency-domain 
received samples and the noiseless samples; and generating a next estimate of the 
starting position of the sequence of desired samples in the received signal sequence 

25 based on an argument of the current frequency-domain correlation between the 

frequency -domain received samples and the noiseless samples. 

In still another aspect of the invention, receiving the sequence of 
desired data symbols further includes using the final improved estimate of the 
starting position of the sequence of desired samples in the received signal sequence 
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to adjust a measurement of frequency offset associated with the sequence of 
desired samples in the received signal sequence. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The objects and advantages of the invention will be understood by 
5 reading the following detailed description in conjunction with the drawings in 

which: 

FIGS, la and lb depict conventional repetition preamble structures; 
FIG. 2a depicts a temporal structure of the discrete-time transmitter 
output samples s[k] during a preamble structure that includes regions G, A and B 
10 in correspondence with those regions shown in FIG. la; 

FIG. 2b illustrates samples received in a receiver after the 
transmitted samples have been transmitted through a noiseless dispersive channel; 

FIG. 3 is a diagram of an exemplary sandamble in accordance with 
one aspect of the invention; 
15 FIG. 4 is a depiction of the temporal structure of a sandamble in 

accordance with the invention; 

FIG. 5a is a diagram depicting hardware components and flow of 
signals and order of operations in a first exemplary embodiment of a first stage 
unit that performs coarse frame synchronization and fine frequency 
20 synchronization in accordance with one aspect of the invention; 

FIG. 5b is a diagram depicting hardware components and flow of 
signals and order of operations in an alternative exemplary embodiment of a first 
stage unit that performs coarse frame synchronization and fine frequency 
synchronization in accordance with one aspect of the invention; 
25 FIG. 6 schematically depicts correlation regions for a received 

signal having G l5 G 2 , G 3 A (including A 2 ) and B (including B 2 ) regions; 
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FIG. 7 is a phasor diagram of a noiseless (coarse) correlation in the 
CF approach and the non-ambiguity region of a fine estimator for the sandamble 
structure; 

FIG. 8 is a graph in which the horizontal axis represents the relative 
5 frequency offsets ^-and the axis is subdivided into the (normalized) ambiguity 

intervals of length D/k 0 ; 

FIG. 9 is a table that presents an overview of the phase angle values 
associated with the frequency offset centroid in the interval A f for various OFDM 
parameter values; 

10 FIG. 10 is a phasor diagram for a relaxed detection with a 

Fine/Discrete technique (F/D) up to | A f | <2, in accordance with one aspect of the 
invention; 

FIG. 11 is a graph in which the theoretical failure rate performance 
of a Coarse/Fine synchronization technique (C/F) is compared to that of the F/D 
15 technique over a range of E s /N 0 values in an Additive White Gaussian Noise 

(AWGN) channel and over a range of relative offset Rvalues; 

FIG. 12 is a graph in which the same failure rate comparison is 
provided for the case in which D c = 10 correlation products are used in the coarse 
correlation; 

20 FIG. 13 schematically depicts the temporal position of the 

cl c2 

associated correlation windows which are used to calculate L [k] and L [k] when 
an alternative sandamble structure having an extended guard interval in a second 
part is used, in accordance with one aspect of the invention; 

FIG. 14 is a graph in which the false-lock rate P ff of CF in an 
25 AWGN channel is shown for various normalized frequency offsets |£ f | e [0.0,0.9]; 

FIG. 15 is a graph depicting the frequency synchronization 
performance of CF for transmission over a severe multipath channel; 

FIG. 16 is a graph depicting the false synchronization probability 
for transmission over the AWGN channel for various values of £ f ; 
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FIG. 17 is a graph comparing simulated versus theoretical results of 
using the FD technique in a multipath channel, assuming transmission over a 
severe multipath channel; 

FIG. 18 is a graph depicting simulated frequency synchronization 
5 failure rates over SNR and theoretical predictions for the inventive FD approach 

with two correlations in an AWGN channel with different a priori probabilities for 
the ambiguity intervals; 

FIG. 19 is a graph depicting the frequency false-lock rate P ff of the 
CF scheme with D c = 6 and various misplacements D m in a multipath channel; 
10 FIG. 20 is a graph depicting the performance of the FD approach 

for equal a priori probabilities of the frequency offset intervals; 

FIG. 21 is a graph depicting the performance characteristic with the 
same a priori parameters previously used; 

FIG. 22 is a graph depicting simulation results for a Fine/Discrete 
1 5 approach with Pilot Correlation (FDPC); 

FIG. 23 is a graph further depicting simulation results for the FDPC 
technique in accordance with the invention; 

FIG. 24 depicts a principle arrangement of differential subcarrier 
pilots within a subcarrier multiplex, in accordance with one aspect of the 
20 invention; 

FIG. 25 depicts an OFDM symbol structure with a prefix and a 
postfix guard interval in accordance with an aspect of the invention; 

FIG. 26 is a diagram depicting hardware components and flow of 
signals and order of operations in an exemplary embodiment of a second stage unit 
25 that performs fine frame synchronization in accordance with another aspect of the 

invention; 

FIGS. 27 and 28 are each graphs showing the simulated standard 
deviation of £ in the iterative frame synchronization in accordance with an aspect 
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of the invention for transmission over AWGN in the presence of an initial 
misplacement £> m ; 

FIGS. 29 and 30 are graphs depicting the time synchronization 
failure rates for transmission over AWGN in the presence of an initial 
5 misplacement Z? m ; 

FIG. 31 is a graph depicting, for postfix symbols comprising 
various numbers of samples, the Packet Error Rate (PER) as a function of SNR for 
8DPSK modulation in the time direction for ideal synchronization and operational 
burst synchronization before and after the frame synchronization refinement; 
10 FIG. 32 is a graph of the PERs plotted as a function of SNR, where 

PERs are determined for the cases of ideal synchronization and operational burst 
synchronization before and after the frame synchronization refinement in 
accordance with the invention; and 

FIG. 33 is a graph showing a comparison between differential 
15 modulation in the frequency direction and in the time direction wherein for each, 

PERs are plotted as a function of SNR for ideal synchronization and operational 
burst synchronization before and after the frame synchronization refinement in 
accordance with the invention. 



DETAILED DESCRIPTION 

20 The various features of the invention will now be described with 

respect to the figures, in which like parts are identified with the same reference 
characters. 

The various embodiments of the invention offer solutions for the 
problem of performing carrier frequency synchronization and frame 
25 synchronization (i.e., with respect to time) of bursty (i.e., spontaneous, single- 

shot) data transmissions received over unknown frequency-selective channels. The 
invention is well-suited for use with radio channels, but may also be applied to any 
other type of channel that exhibits similar characteristics. In the following 
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discussion, an overview of the invention is first presented. This is followed by a 
more detailed description of the various aspects of the invention. 

In one aspect of the invention, the synchronization procedure 
comprises two stages, with the first stage itself comprising two steps. In the first 
5 step of the first stage, a coarse timing estimate and a (low-variance) frequency 

offset estimate are extracted. The latter can be corrected even though a frequency 
ambiguity may occur. The second step in the first stage aims to resolve this 
frequency ambiguity so that the actual carrier frequency offset can be estimated 
over a wide range. 

10 The second synchronization stage is finally dedicated to improve the 

initial timing estimate. Depending on the specific synchronization problem setting, 
each of the two stages can be used in a stand-alone mode or, alternatively, they can 
be applied jointly. 

The first stage performs a time-domain processing (e.g., 

1 5 correlation) of samples to exploit a periodic signal repetition and to extract the 

coarse timing, frequency offset and to resolve frequency ambiguities. Here, 
another aspect of the invention divides a preamble portion of the transmitted signal 
into two widely spaced signal parts. For example, a data signal (or other type of 
signal) may be embedded in-between the " preamble" parts to form what is referred 

20 to herein as a Sandwich Preamble (SP), also referred to herein as a "sandamble". 

The use of the inventive sandamble is very advantageous in terms of frequency 
accuracy. 

The second stage estimates the fine time offset of the received 
modulation signal. A frequency-domain approach is chosen that operates on some 
25 frame of signal samples that has been identified as the Orthogonal Frequency 

Division Multiplexing (OFDM) symbol of interest by some coarse time estimator. 
The initial timing estimate may stem from the first synchronization stage or it may 
alternatively be the result from relatively fixed periodic communication protocol 
structures in the specific transmission system. For this second synchronization 
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stage, the synchronization preamble is required to be an OFDM symbol. From the 
spectrum of this signal part, an improved time offset estimate can be obtained by a 
correlation technique across the subcarrier amplitudes. If the OFDM parameters 
(e.g., number of carriers, type of subcarrier modulation) permit, a decision- 
directed version of the algorithm can even estimate the time offset from any 
information-carrying OFDM symbol in the burst, so that no additional redundancy 
is required to estimate the offset up to a maximum value given by the specific 
parameter set. That is, the modulation scheme for data transmission in the burst 
need not generally be restricted to the use of OFDM, but it is advantageous to do 
so to permit the second synchronization stage to use a decision-directed algorithm 
to estimate the time offset from an information-carrying symbol in the burst. 

Consequently, data-aided as well as decision-directed estimation can 
be used in the second stage, wherein data-aided estimation requires additional 
redundancy but offers more robustness. Based on the steadily improved estimate 
of the time offset, the time window may be iteratively realigned until the estimate 
converges to an optimum estimate with minimum variance. The coarse initial time 
offset estimate need only be reliable enough to ensure convergence. 

For operation over a multipath (ISI) channel, the guard interval is 
advantageously split into prefix and postfix parts, as the time synchronization from 
the proposed estimator is systematically delayed by the center of power gravity 
time of the respective power delay channel profile. Hence, the estimated 
synchronization positions of the estimator are usually shifted to positive values. 

Having presented an overview of various aspects of the invention, a 
more detailed description of these and other aspects will now be presented. 

Because the invention uses a new type of repetition preamble to 
perform both time (frame) and frequency synchronization of a bursty OFDM 
transmission over unknown frequency-selective channels, this will first be 
discussed. 
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Signal Repetition 

As mentioned in the BACKGROUND section, a repetition preamble 
is commonly proposed as a burst training sequence in order to allow frame- and 
carrier frequency synchronization in digital transmission over unknown and 
5 severely dispersive channels. The general principle of exploiting a cyclic training 

signal for frame synchronization and carrier frequency offset estimation was 
originally suggested for single-carrier transmission in [CMU87]. Linear as well as 
memoryless non-linear channel distortions will have negligible effect on the 
synchronization performance. Thus, the core of the investigated coarse frame and 
10 carrier frequency algorithm is to exploit two regions of signal repetition, each 

consisting of a useful synchronization length of modulation intervals. The 
two regions are spaced k 0 samples apart. The parameter k 0 is denoted as the 
correlation basis. 

There are several possibilities of how the OFDM transmitter can 
15 introduce the signal repetition: 

Extended guard interval: This possibility applies for very large 
OFDM symbols (e.g., Z)=512), so that the "natural" periodicity 
due to the guard interval is slightly extended by some (e.g., 20) 
additional guard samples to allow for synchronization. In this 
20 example, we would have a synchronization length of D^ c =20 and 

a correlation basis of k 0 =D— 512. 

OFDM symbol repetition: This possibility is preferable for small 
OFDM symbols (D^128) and this means that an entire OFDM 
symbol of dimension is repeated to allow for synchronization. 
25 The value for the correlation basis k 0 can be varied. 

The OFDM symbol that is repeated to form the preamble is 
used to carry data, too. Hence, the preamble structure itself is - 
apart from the periodicity - random and only one half of the 
preamble samples represent training overhead, even though this 
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additionally provides a time diversity factor of 2 for data 
demodulation. 



10 



Mathematical Description 

For the mathematical notation we assume a non-oversampled 
OFDM transmitter and receiver. Furthermore, we assume a perfect match of 
sampling frequencies in both the transmitter and the receiver, resulting in 

r T =r R =r. 

The exploitable time-domain property in the transmitted 
synchronization structure is the repetition property 

s[k + k 0 ] = s[k] 



(1) 



Note that without restriction on generality, the beginning of the first 
useful synchronization symbol was chosen to be k=0. For negative k, the guard 
interval of the synchronization sequence is transmitted. 

For the continuous-time transmit signal, one obtains 



s(t) = 



(2) 



15 



After being subjected to a dispersive channel, additive noise, a 
receive filter characterized by h R (t), and a carrier frequency offset is 

t \ 



(3) 



with h{t) being the overall impulse response defined as h(t) = h 7 (t)*h c (t)*h R (t) . 
where h^t) is the impulse response of the transmit filter and is the impulse 
response of the channel, and "*" is the convolution operator. 
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The noiseless received signal clearly is 



fj(0=^ A/ "'-£ s[klKt-k'T). (4) 



The noiseless signal is sampled so as to yield 



where we assume that | Ar so | «r and that the number D e of sample indices k, where 
h(kT) is significantly non-zero does not exceed the guard interval of length D g . 
Hence, we obtain a received sample sequence exhibiting the signal property 



r[k+k 0 ]=e v2 ^ T -f[k] Vk e[-(D g -D c ), D^-l] (6) 

at least in the mentioned interval. This is true if the guard interval is not shorter 
than the channel impulse response, i.e., Dg>D Q is valid. Clearly, the phase 
relation between £ 0 -spaced samples is independent of Ar so . 

The systematic phase increment of 2nAf co Tk 0 over an interval of k 0 
samples expressed in Eq. (6) can be exploited to estimate the frequency error at the 
receiver site. This will be discussed in further detail later in this disclosure. 

With the OFDM subcarrier spacing A/ su5 = l/(DT) and the 
Normalized Carrier Frequency Offset (NCFO) £/=(A/ co )/(A/ sub ) we may write 
27cA/ co r=(27i/Z))^ and hence we can rewrite Eq. (6) to yield 

r[k+k 0 ]=e VD °-r[k] V* e[-(Z> g -Z> e ), D^-l] ( 7 > 
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Properties of the Received Signal 

Apart from a phase rotation, the received noiseless signal exhibits a 
repetition with the periodicity interval k 0 within a region of at least D sync samples. 
Thus, a metric that reflects the similarity of the signal spaced k 0 samples apart can 
5 be applied in the receiver to reliably detect the presence of a synchronization 

sequence. Suitable metrics are discussed later in this disclosure. 

If D G <D g a portion of D g -D e samples of the guard interval will be 
unconsumed. (If h c (t) has only small magnitude values towards the end of the 
impulse response, a corresponding portion of samples will not be identical, but 

10 will nevertheless be "very similar" so that the "unconsumed" part seems to be 

even longer.) The awareness of this property is important for understanding the 
arising problem of uncertainty in the time synchronization. The position of the 
symbol start cannot be detected with absolute exactness by the maximum metric 
detection proposed in [SC96]. Hence, only a coarse estimate for the timing offset 

15 is obtained. 

The robustness for detecting the presence of the synchronization 
structure and the accuracy of the frequency estimate is dependent on the parameter 
D S y nc because this one determines a "processing gain" and thus the noise 
resistance of the approach. Furthermore, fc 0 influences the frequency estimation 
20 accuracy and the permissible frequency range. By variation of D S y nc and £q, the 

synchronization approach with signal repetition becomes scalable in the achievable 
accuracy. 

Requirements for Applying this Synchronization Structure 

Clearly, the duration of the exploited synchronization sequence may 
25 not exceed the coherence time of the channel. This is usually not a critical 

requirement in intended applications. 
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It must be provided that the maximally occurring frequency offset 
does not exceed the number of virtual carriers on either side of the transmit 
bandwidth. In other words, the received signal must "hit" with all active 
subcarriers into the passband of the receive filter H R (f) without prior frequency 
5 control. Thus, all desirable receive signal information is contained in the 

Equivalent Complex Baseband (ECB) domain. After the frequency offset has been 
sufficiently corrected, it is possible to demodulate the frequency multiplex via a 
Discrete Fourier Transform (DFT) without suffering from excessive interference 
from other subcarriers; the frequency correction must approximately restore the 
1 0 mutual subcarrier orthogonality . 

Conventional Repetition Preamble 

For these examples, it will be assumed that the number of carriers 

in the system is relatively small, so that the discussion will focus on the repetition 

of a full OFDM symbol of dimension Z> sync . Thus it will be assumed, for the 

15 purpose of example, that the preamble comprises one repeated data-carrying (and 

therefore random) OFDM symbol. (Of course, the preamble may alternatively 

comprise known data.) Several conventional repetition preambles were discussed 

in the BACKGROUND section of this disclosure, and illustrated in FIGS, la and 

lb. Examining the conventional preamble structure of, say, FIG. la more closely, 

(2) 

20 it can be seen that k 0 =D sync and that £> g v =0. We arrive at two identical and 

directly successive OFDM signal parts that we call the P structure. Such a 
synchronization structure consists of 2D sync time-domain samples plus D g samples 
in a guard interval. As usual, the guard interval consists of a cyclic prefix so that 
one small part of the transmit sequence occurs three times in the synchronization 

25 sequence. The D g last samples of the first D sync -carrier synchronization OFDM 

symbol simultaneously represent the "guard" interval for the second OFDM 

symbol because both subsequent symbols are identical. That is the reason why we 

(2) 

were allowed to set D g —0. This decisive property of the P structure is depicted 
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schematically in FIG. 2a. In particular FIG. 2a depicts a temporal structure of the 
discrete-time transmitter output s[k] (each transmitted sample being a complex 
value) during the preamble structure P that includes regions G, A and B in 
correspondence with those regions shown in FIG. la. In each region of this and 
5 all subsequent figures in this disclosure, the triangular shape is chosen to illustrate 

the sequence of samples in one region (e.g., the A region) and its repetition in 
other regions (either complete repetition as in the B region's repetition of the 
samples first transmitted in the A region, or partial repetition as in the G region). 
The triangular shape is not intended to mean that the power of the samples 

10 decreases towards the end. For example, by the triangular shape, one can see that 

the set of samples in the B region are identical to those transmitted in the A region, 
and that the set of samples in the G region are identical to a rightmost part of each 
of the A and B regions, thereby providing cyclic repetition. 

FIG. 2b illustrates the received samples, r[k], in the receiver after 

15 the transmitted samples, s[k] have been transmitted through a noiseless dispersive 

channel with either linear or nonlinear distortion and frequency offset Af co . 
Synchronization is based on these received samples r[k]. It can be seen that the 
received guard interval, G\ is not entirely consumed by the channel. With this 
conventional repetition preamble, choosing D sync =D/2 is favorable from a 

20 hardware design point of view because this permits the OFDM symbol structure of 

the transmission system to be the same for all symbols. 

Sandwich Preamble (Sandamble) Structure 

In accordance with one aspect of the invention, a novel type of 
repetition preamble structure, called a Sandwich Preamble (henceforth, 
25 "Sandamble"), is utilized. An exemplary sandamble is depicted in FIG. 3. The 

exploitable periodicity intervals are illustrated by the lines 301, 303 and 305. The 
synchronization structure has the following form: two identical and repeated signal 
parts that are separated from one another by one or more regular OFDM symbols 
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of the burst. In the exemplary embodiment of FIG. 3, the set of samples contained 
in the B region are identical to those contained in the A region. Moreover, three 
guard regions are provided, here designated (Gj, G 2 ) and G 3 . The notation in 
which parentheses enclose the designations of Gj and G 2 is used to indicate that 
5 the first and second guard regions together comprise an extended guard region. 

The guard regions G± and G 2 are contiguous to one another, and together precede, 
and are adjacent to, the A region. The third guard region, G 3 , precedes and is 
adjacent to the B region. The third guard region G 3 is separated from the A region 
by one or more regular OFDM symbols of the burst. 

10 The signal in each of the extended and third guard regions, (G lf G 2 ) 

and G 3 is a copy of a rightmost part of the signal included in the A region (i.e., 
the part of region A designated A2 in FIG. 3). (Of course, the signal in each of 
the guard regions (G x , G 2 ) and G 3 may also be considered a copy of the rightmost 
part of the signal included in the B region, since the signal in the B region is, 

15 itself, a complete copy of the signal in the A region.) Since the length of the 

extended guard region, (G lf G 2 ), need not be the same as the length of the third 
guard region, G 3 , different length copies of rightmost parts of the signal in the A 
region are used. 

In the sandamble, the first guard region G 1 serves the same purpose 
20 as the guard region G found in the conventional repetition preamble depicted in 

FIG. la. The second guard region G 2 serves to extend the guard region beyond 
G^ The third guard region G 3 is additionally necessary for use as a prefix for the 
repeated part B due to the split-up (i.e., separation) of the preamble components. 
FIG. 4 is another depiction of the temporal structure of a 
25 sandamble. Here, the triangular shape is useful for depicting which portions of the 

signal are replicated in various parts of the sandamble. It can be seen from FIG. 4 
that the burst frame starts with the synchronization OFDM symbol having a sample 
length of £> sync , preceded by an extended guard interval of length D^ l \ Then, 
regular (i.e., information-carrying) OFDM symbols with D carriers together with 
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regular guard intervals of Z> g samples may follow. Finally, the second 
synchronization OFDM symbol (identical to the first) of length D sync together with 
a guard interval of length D^ 2 \ not necessarily equal to D^ l \ marks the end of 
the sample frame (partly) used for synchronization parameter estimation. 
5 Thereafter, the remaining regular (information-carrying) OFDM symbols follow in 

the burst. 

By splitting up and enclosing regular OFDM symbols, the 
correlation basis £ 0 is enlarged. This quadratically reduces the estimation variance 
of the carrier frequency estimate, while also linearly reducing the 16ck-in range. A 
10 specified maximum tolerable estimation variance can be achieved with significantly 

reduced synchronization overhead (dedicated training samples), when k 0 is large. 

In the following, the sandamble is presumed to have just one 
embedded regular OFDM symbol. 

Conventional Preamble vs. New Sandamble 
15 A first disadvantage of the sandamble is the reduced frequency 

carrier offset lock-in range. However, in practice the frequency carrier lock-in 
range usually does not need to be very large. Thus, in accordance with another 
aspect of the invention, two methods to circumvent ambiguity by a second 
operation step on a slightly enlarged guard interval of the first 
20 synchronization OFDM symbol will be described later in this description. For two 

of the presented frequency offset ambiguity resolution algorithms, D^>D g is 
required. 



A second disadvantage of the sandamble is that in comparison to the 

conventional repetition preamble, the second part of the synchronization preamble 

25 requires an own guard interval, so that the total "preamble length" is slightly 

increased from D g + 2D sync in the conventional repetition preamble to 

D rt ^+D„^+2D Qvnr in the sandamble. If it is assumed that the usual choice will 

g g. by lit- 

he dJ®~D „, then the number of additionally required modulation intervals for 
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synchronization is Dg l \ The sandamble structure suggests greater efficiency than 
the conventional repetition preamble when the achievable frequency estimation 
variance o^ f is compared with the transmission overhead. That is, the sandamble 
structure achieves the best accuracy (low variance) if the number of training 
samples (overhead) is fixed. 

FTRST STAGE: COARSE TIME AND FINE FREQUENCY 
SYNCHRONIZATION 

As mentioned earlier, the synchronization procedure in accordance 
with the invention comprises two stages, with the first stage itself comprising two 
steps. In the first step of the first stage, a coarse timing estimate and a low- 
variance frequency offset estimate are extracted. The latter can be corrected even 
though a frequency ambiguity may occur. The second step in the first stage aims 
to resolve this frequency ambiguity so that the actual carrier frequency offset can 
be estimated over a wide range. 

A block diagram of a first exemplary embodiment of the first stage 
is depicted in FIG. 5a. This figure depicts not only the hardware components used 
in the exemplary embodiment, but also the flow of signals and order of operations, 
and as such also serves as a flowchart of the various steps performed in the 
embodiment. In order to depict this flow of operation, the same hardware 
component, such as the sample memory 501, may in some cases be depicted in 
several places in the figure. Notwithstanding this multiple depiction, only one of 
any such unit is required. The functions of the various units described in this and 
subsequent figures may be implemented in any of a number of ways. For 
example, programmable processing logic may perform the herein-described 
functions by executing program instructions that have been stored in a suitable 
computer-readable storage medium (e.g., Random Access Memory (RAM), 
magnetic storage medium, optical storage medium, and the like). Alternatively, 
hardwired logic may be designed and built to perform some or all of the various 



WO 99/65180 



PCT/EP99/03793 



-28- 

functions. Each of these alternatives, as well as their equivalents, are intended to 
be considered within the scope of the invention. It will further be understood that 
the segregation of the various functions into the illustrated blocks is for the 
purpose of facilitating a description of the invention. In practice, it is possible for 
5 some or all of the illustrated blocks to be combined into a single unit that performs 

all of the combined functions. 

Turning now to a discussion of FIG. 5a, this first stage does not 
restrict the type of modulation used in the preamble symbol. Single-carrier 
modulation as well as OFDM symbols can be used as a preamble. 

10 Received (noisy) samples r[k] are supplied to a sample memory 501. 

It is known that the stored received samples r[k] will include a sandamble structure 
503. In a first step of the first stage, the received samples r[k] are supplied to a 
first correlation unit 505 that generates a coarse frame start (time) estimate 509 as 
well as a first frequency correction 507 that is applied to the received samples r[k]. 

15 These are generated by determining similarity metrics from the received samples 

r[fc], and using information obtained from the generated metrics to detect the 
sandamble 503. From the detected sandamble 503, both the first frequency 
correction 507 and the coarse frame start (time) estimate 509 may be determined. 
Alternative embodiments of the invention may be derived by using different types 

20 of similarity metrics. Suitable metrics will now be described. 

Minimum-Mean-Squared-Error (MMSE) Criterion 
To approach the problem of frame and frequency synchronization, a 
sequence of D sync contiguous noisy received samples is collected in the vector 



r k =(r[k],...s[k+D-\]) r 



(8) 
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The (fictive, and therefore non-observable) noiseless received 
samples are collected in 



(9) 



Finally, with the definition of the noise vector 



¥ («M v ..,«[^v-i]y 



(10) 



we have the relation r k = r k +n k - 



The synchronization is based on maximizing the similarity 
probability of sample sequences. A valid time and frequency synchronization is 
achieved if the synchronization tuple (fc, ^ is contained in the valid 
synchronization tuple set (valid region) 



~ 



(id 



All time positions that exhibit the periodicity property, i.e., the 



identity 



v -2jt— c; 

between r. Ulc and e D f k , are valid. In the presence of an 



unconsumed guard interval the valid k' will form contiguous intervals. Apart from 
that, the NCFO is a valid estimate for all ^.=^+Xy- with x e Z , i.e., there 
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is an NCFO ambiguity interval of D/k Q . An unambiguous NCFO estimate is a 
priori impossible if the frequency parameter to be estimated is not restricted in 
range. Consequently, there may exist a multiplicity of valid synchronization 
tuples, so that the cardinality of R k ^ is larger than one. 

For periodicity detection we introduce the decision vector 



d- = r- -e**&r- < 12 > 



V2^lj ^4«/ (13) 



where for the first two terms in Eq. (13), it is noted that 

r^ Q - e V2 ^ f r f = 0 V^e^. 



The decision vector given by Eqs. (12) and (13) can be used to test 

the two synchronization hypotheses £ and £ simultaneously. It provides 

2 7? 

reasonable information at least for o s > a n . In the case of valid hypotheses, 

d^ represents a complex- valued zero-mean D^^-dimensional Gaussian 
distributed random variable. This property follows directly from Eq. (13). In the 
noiseless case d £ % will be exactly zero for the ideal frame positions and the 
perfect NCFO estimates . In the presence of mutually uncorrelated white Gaussian 
noise vectors n £ ^ and the noise in d^ will again be white Gaussian with 
a noise variance of 2o n per complex dimension. 
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The probability density function (pdf) of d^ under the condition 
of a valid synchronization tuple is 



fa) ^ *J = 



II d. 



\\2\ 



exp 



2ol 



(14) 



where || d \\ 2 =d H d . Clearly, (•) denotes the complex conjugate transpose of a 
vector. 

The frame offset and the NCFO have to be estimated jointly. The 
joint estimates are obtained by performing 




The argmax operator yields the argument (tuple) that maximizes the 
given expression. 

This joint synchronization estimate is equivalent to 

£JL = argmin \\d r A\ 2 = argmin \\r~ - e D r&\ 2 . (16) 
V 1 U,) ^ (M f ) 



The argmin -operator yields the argument (tuple) that minimizes the 
given expression. The expression to be minimized can be modified to yield 



11%/ = \h, k f + \htf ~ 2St\e D r»r £ ^ 



= P[k+k 0 ]+P[k] - 23t[e D S[k]\ , 



(17) 
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where we introduced the complex correlation 

D -1 



S[k] = £ r*[£+K] • r[k+k Q +K] 



(18) 

K=0 



and the power sum inside a frame of subsequently received samples 

P[k] = J |r[*+K]| 2 . (19) 



Both values, S[k] and P[£], can be recursively calculated from 
S[k - 1] and P[k - 1], respectively. This is accomplished by adding one new 
element and subtracting the oldest one. Under the assumption of negligible 
rounding errors, the sum need not be calculated from scratch in each modulation 
interval. 

The expression in Eq. (16) is exactly the same metric which was 
proposed for use as a periodicity metric for joint frame- and frequency 
synchronization in [CMU87]. It is a norm for the Mean-Squared Error (MSE) 
between the received samples, spaced k 0 samples apart and can therefore be 
exploited to monitor the degree of periodicity in the signal. The joint timing and 
frequency offset decision in the two-dimensional metric in Eq. (16) is in favor of 
the frame start hypothesis k and the NCFO hypothesis £ j which achieve the 
minimum MSE (MMSE). 

With the modified metric 

M[k] ^ P[k+k 0 ] + P[k] - 2\S[k]\ (20) 
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the original two-dimensional search can be broken down into two one-dimensional 
estimation problems. The estimation of the frame position comes first via the 
minimum MSE criterion 

k = argmin M[k] m 



Thereafter the Maximum Likelihood (ML) estimate for the NCFO 
can be obtained by evaluating S[k] at the estimated frame start, i.e., k = k, and 
this yields 

^ = arg ^ ■ (22) 



Clearly, |^| <D/(2k 0 ) is a minimum requirement for the 
unambiguousness of the frequency offset estimate ^ The quality of this estimate 
is further investigated later in this description, but it is mentioned here that for ISI- 
free frequency offset estimation, £e[-(D g - £> e ), 0] must be ensured. 

Maximum-Likelihood Criterion 

In [SvdBB95] and [vdBSIB95], Sandell, van de Beek and Borjesson 
propose an optimum metric that is based on a ML frame synchronization approach. 
The detailed derivation can be found in Jan-Jaap van de Beek, Magnus Sandell, 
and Per Ola Borjesson, ML Estimation of Timing and Frequency Offset in 
Multicarrier Systems , Research report, Div. of Signal Processing, Lulea 
University of Technology, Sweden, 1996 (henceforth, "[vdBSB96]"). The 
received signal is modeled as a complex zero-mean Gaussian distributed white 
random process, which is only valid for some OFDM transmit signals. Unlike the 
case with large numbers of unused subcarriers , as described in Stefan Muller and 
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Johannes Huber, "A Novel Peak Power Reduction Scheme for OFDM", 
Proceedings of the International Symposium on Personal. Indoor and Mobile Radio 
Communications (PIMRC97) . pages 1090-1094, Helsinki, Finland, September 
1997 (henceforth, "[MH97b]"), or with adaptive modulation, the zero-correlation 
5 of samples is not generally true for OFDM transmit signals with strongly varying 

transmit powers in the subcarriers. It is especially not true for the received signal 
because it was convolved with some dispersive channel impulse response. 
Consequently, the derivation in [vdBSB96] is itself based on the assumption of a 
non-dispersive channel. In the course of their analysis, they first arrive at a two- 
10 dimensional metric for joint time and frequency estimation and with the same 

reasoning as in the previous subsection they derive the one-dimensional search 
criterion 

k = argmin (p(P[£+* 0 ] + P[k]) - 2 \s[k]\j (23) 



for frame synchronization, where the constant 

2 



P = ~ 2 (24) 



accounts for the Signal-to-Noise Ratio (SNR) at the receiver input. Apart from 
15 this SNR-adaptive factor the metric has the same structure as the one in Eq. (20). 

Clearly, P[k] and S[k] are defined in Eqs. (18) and (19). Even though the 
derivation assumes a non-dispersive channel, the metric retains its optimally and 
superiority in dispersive channels. 



20 



Maximum-Correlation Criterion 

A simplified frame synchronization metric is proposed in T. Keller 
and L. Hanzo, "Orthogonal Frequency Division Multiplex Synchronisation 
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Techniques for Wireless Local Area Networks", Proceedings of the International 
Symposium on Personal. Indoor and Mobile Radio Communications (PIMRC96V 
pages 963-967, Taipei, Taiwan, 1996 (henceforth U [KH96] W ). Here, we obtain the 
frame start via 



which represents the time position of Maximum Correlation (MC) magnitude (cf . 



square roots need to be processed in an implementation. 

If the received signal has a constant envelope and the noise is 
moderate so that P[k] = const V£, then the latter criterion would be as optimal as 
the criteria in Eqs. (23) and (21). However, this is definitely never true for 
OFDM signals because the signal envelope of OFDM is far from constant. It is 
not even constant for multipath-corrupted receive signals in single-carrier 
modulated systems. Consequently, the criterion in Eq. (25) must be suboptimal, 
because it does not account for the average power inside the currently processed 
synchronization window. 



k = argmax \s[k]\ , 



(25) 



Eq. (18)). Equivalently, the maximum of \s[k]f 



can be the criterion, so that no 



A Fourth Criterion 



In [SC96], Schmidl and Cox suggest applying the "defined 1 [SC97] 



metric 



k = argmax 



(27) 




for which the required values are defined in Eqs. (18) and (19). 
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A Fifth Criterion: Maximum Normalized Correlation 
The inventors of the herein-described and claimed invention have 
developed yet another criterion, referred to herein as Maximum Normalized 
Correlation (MNC) that, they consider to be preferable for use in performing coarse 
frame synchronization with the new sandamble. The metric in accordance with 
this aspect of the invention, and its use in determining a coarse time 
synchronization estimate (i.e., k) is 



k = argmax [ L Jl — • ^8) 



k 



(P[k+k Q ]+P[k]f 



In operation, then, the first correlation unit 505 (see FIG. 5) 
determines, for each of a number of hypothesized coarse frequency estimates, £, a 
10 correlation value representing the correlation between the hypothesized A 

region and the hypothesized B region of the received samples r[k] . These 
correlation values may be stored in a correlation memory 511. After computing 
all of the correlation values, the first correlation unit 505 then examines these to 
perform the argmax-operation, that is, to find that value of k that maximizes the 

15 metric . This maximizing value, £, is the coarse frame start 

(P[k+k 0 ]+P[k]Y 

estimate 509. 



Reducing Frequency Estimation Variance 

Having estimated the locations of the A and B regions in the 

received samples r[k], it is further possible for the first correlation unit 505 to 

20 estimate the frequency of the received signal. The prevalent problem in frequency 

estimation is the frequency accuracy. The variance of the frequency estimate is 

2 

inversely proportional to k 0 . A second problem is the estimation range which 
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gives the maximum allowable frequency offset which can be recognized by the 
estimator without ambiguity. For each of the preamble structures depicted in 
FIGS, la, lb and 3, the estimation range for the frequency offset derived from the 
correlation between the A and B regions is inversely proportional to £q. It 
5 follows, then, that an enlarged periodicity interval has the effect of reducing the 

estimation range. 

A number of possibilities can be outlined to obtain accurate (i.e., 
low-variance) frequency estimates. These are: 

With reference to the repetition preamble depicted in FIG. la, an 
10 increase in the number of samples within the repeated signal parts in 

the A and B regions improves the estimation accuracy. See [SC96]. 
The effect is that the diversity factor (and the training overhead) 
becomes larger and the periodicity interval is enlarged. Hence, the 
estimation range is reduced. 

Using several repetitions instead of only two gives the opportunity 
to use various periodicity intervals, as is depicted in the preamble 
structure of FIG. lb. Here, the processing of the regions A and B 
as well as the regions B and C permits larger frequency offsets to be 
estimated, while the processing of the regions A and C results in a 
better frequency variance. See [FM93] and [FM94]. 

In accordance with one aspect of the invention, another approach to 
obtain accurate frequency estimates is to increase the periodicity 
interval without increasing the number of repeated samples. Here 
the diversity factor (and the training overhead) remain 
25 approximately the same, while the estimation range is similarly 

reduced. By doubling the periodicity interval k 0 between the A and 



15 



20 
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B regions, the estimation variance is reduced by a factor of 4, while 
the estimation range is reduced by a factor of 2. This idea of 
trading variance against estimation range is realized via the very 
flexible sandamble approach, in which the repeated synchronization 
symbols are separated by data. 

In terms of the relationship between frequency estimation variance 
and training overhead, the sandamble has considerable advantages over the pure 
(conventional) preamble types. This is already apparent in FIGS, la, lb and 3 by 
comparing the number of training samples. Depending on the specific parameters, 
the number of overall training samples in the sandamble can be half that of the 
preamble types with comparable frequency estimation variance, but clearly with a 
reduced estimation range. The latter can again more cheaply (in the sense of the 
number of dedicated samples) be enlarged by applying the method in the next 
Subsection. 

Resolving Frequency Ambiguities 

This discussion presents an overview of the approach to resolving 
frequency ambiguities that is used by the first stage of the exemplary embodiment. 
Following this overview, the various aspects will be examined in greater detail. 

Returning now to the sandamble depicted in FIG. 3, the extended 
guard region G 2 is used to resolve the frequency ambiguities which may happen 
when correlating regions A and B. In particular, the first frequency correction 507 
generated from the correlation of the A and B regions by the first correlation unit 
505 is used to correct (step 513) the frequency offset of the received samples r[k]. 
The frequency correction step 513 may be performed, for example, by multiplying 
samples (stored in the sample memory 501) by a complex-valued rotation factor 
that is a function of the estimated frequency offset value (e.g., e -j2n&jrk, where &f 
is the estimated frequency offset). This frequency offset may be inaccurate, 
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however, because of the possible occurrence of a periodic frequency ambiguity 
with a parameter-dependent interval. Consequently, this residual discrete- valued 
frequency offset needs to be detected. It is important to realize that there is only a 
discrete set of valid residual frequency offsets, so that this procedure is a detection 
rather than an estimation. 

In order to obtain a decision variable, G 2 is correlated with A 2 by a 
second correlation unit 515. In addition, a third correlation unit 517 correlates the 
G 2 region with the B 2 region. In one aspect of the invention, the first frequency 
correction 513 makes it possible for the two correlation results (from the second 
and third correlation units 515, 517) to be combined (e.g., by addition) by a 
combiner 519, although such combination is not essential to the invention. That 
is, in alternative embodiments, only one or the other of the two mentioned 
correlations need be performed. Returning now to the exemplary embodiment in 
which both correlation results are generated and combined, it is found that this 
combination yields an additional 1.55 dB gain for the frequency interval decision. 
From this combined correlation result, the residual carrier frequency offset is 
detected (not estimated) by a frequency offset interval detection unit 521. This 
frequency interval decision 523 is combined with the first frequency correction 
estimate 507 from the first correlation unit 505 to form an overall frequency offset 
estimate 525 (e.g., by multiplying the stored sample by a complex-valued rotation 
factor that is a function of the estimated frequency offset value). Because this 
aspect of the invention combines the making of a fine estimate with the making of 
a discrete decision, this frequency offset estimation procedure is herein referred to 
as "Fine/Discrete" (FD). 

A second advantage of the FD approach lies in the non-uniform 
distribution of frequency offsets in practical transmission systems. Usually, the 
frequency offsets around zero are more probable than the large carrier frequency 
offsets. In this case, the FD approach allows the use of a priori probabilities of 
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frequency offset intervals in the discrete interval decision. This minimizes the 
overall frequency false-lock probability of FD. 

Having described an overview of the operation of the first stage of 
the synchronization procedure, the various aspects of this first stage will now be 
5 described in greater detail. 

Three Approaches for Resolving Frequency Ambiguity 

In this section three approaches for avoiding or resolving frequency 
estimation ambiguities are analyzed. This may increasingly occur with the new 
sandamble structure if the frequency offsets are higher than or near the lock-in 

10 range limit of the sandamble. In the first two of the following three techniques, 

one of the two guard intervals of the two synchronization symbol parts is 
extended, so that correlations with different basis lengths become feasible. By 
appropriate processing of the correlation results, the ambiguity can be resolved or 
avoided in certain limits, thereby leading to a rather feed-forward kind of 

15 frequency estimator. Referring back to the preamble structure of FIG. 4, the 

reader is reminded that in the exemplary sandamble, >D^ =D^. At this 
point the variable D c is further introduced, which is the number of correlation 
products that, due to the guard interval extension, can be used to resolve frequency 
ambiguity. Clearly, D C <,D^ - D e must be provided in order to ensure the 

20 required periodicity for the correlation products. 

The third approach does not require a guard interval extension and 
is therefore more efficient in throughput. It exploits the knowledge of differential 
subcarrier pilot symbols, introduced for the (iterative) frequency-domain time 
synchronization. In order to determine the frequency offset, a trial-and-error 

25 technique is used, which is computationally more expensive than the first two 

approaches. 
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First Technique: Coarse/Fine (CF) Approach 
In this technique, ambiguity is avoided prior to the use of the 
sandamble approach for performing the fine frequency offset estimation. A similar 
synchronization scheme is described in [LSM97] . In a preprocessing unit, a 
coarse frequency estimate is obtained from a first correlation between the A and B 
regions of the preamble structure shown in FIG. lb. This first correlation has a 
shorter correlation basis. The hereby estimated (potentially large) frequency offset 
is corrected prior to the fine estimation according to the result from a second 
correlation obtained between the A and C regions of the preamble structure shown 
in FIG. lb. 

In accordance with one aspect of the invention, this technique may 
be applied to the new sandamble structure, which permits both the first and second 
correlations to be performed between components having respective shorter and 
longer correlation bases. 

The correlation to use for coarse frequency offset estimation 
exploits a basis of £> sync , consists of D c products and is defined by 

L c [k] = £r*[*-KMA;+Z^ c -K] . (29) 



In FIG. 6, this correlation is schematically depicted for a received 
signal 601 having G!,G 2 , G 3 A (including A 2 ) and B (including B 2 ) regions. 

By substituting, in place of the noisy sample, the sum of the 
noiseless sample plus noise, i.e., r[k] = r[k] +n[k] , we approximately yield 

L c [k] * £ n*-K]rl*+D^ c -K] 

K ~\ (30) 
+ Er[M«[^ c -K] + n'ik-mk+D^-K] , 
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where the noise-noise product has been immediately neglected. Now, we assume 
perfect time synchronization so that L°[k] is evaluated at the correct time instant 
(*=0) and yield 



L c [0] * £|r[-K]| 2 e +J2 * Af " D »™ T 

K = l 

* En-K] 



r*[-K] 



(31) 



10 



where we exploited the property r[D^ nc -K]=f[-K]e +J2nAf ^^ T , 1<;k<D c for the 
calculation of the useful part of the correlation result. Note that due to the same 
property the multiplicative factor for the noise component n *[-k] , namely the 
ratio of the noiseless received signal samples r[D -k] over r*[-K] , has 
unity magnitude. 

At this point, we introduce the new random variable u which 
represents the sum of useful power in the exploited signal part. It is defined as 



U = £ft-K]| 2 , 



(32) 



so that the useful component L c [0] of L c [0] in Eq. (31) can be expressed as 
L c [0] =u-e +J2nAf ^" cT . The variance of the noise component in L c [0] is 



-/T[-k] 



D 



K = l 



f*[-K] 



E kl-K]| 2 -h[^ c -K] + - J n *[-k] 



r*[-K] 



(33) 



K^„ c -^3 + ^ ^" c . j «'[-k] 
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where we assumed that the received signal samples and the noise samples as well 
as the single noise samples are statistically independent. Note that the complex- 
valued noise component in L c [0] is again Gaussian distributed with variance 

o 2 
w-2o n . 

Now, we focus on the statistical properties of the random variable u 

2 

in Eq. (32). If D c ~ D sync , we will constantly find that u = D sync o s , which is 
ensured by the theorem of Parseval. Thus, u is not a random variable for the case 
D c = D S y nc . But owing to preamble efficiency, the more prevalent case 'will be 
D c < <D sync and this will be assumed in the following derivation. It follows 
directly from the properties of the OFDM transmit signal that the noiseless time- 
domain samples r[k ] can be approximately interpreted as complex-valued 

uncorrelated and Gaussian distributed (statistically independent) random variables 

2 

with zero mean and variance given by o s . Hence u is centrally chi-square 
distributed with 2D C degrees of freedom. The pdf of u is then 

\ 



PS U ) = —r- D u exp 



u 



S_j(w) , (34) 



where 6 j (w) is a unit step function that is equal to 1 for all u greater than or equal 
to 0, and equal to 0 otherwise. 

The following analysis aims to find an analytical expression for the 
probability of frequency estimation failure, i.e., that the coarse estimator is not 
able to provide a coarse estimate such that the fine estimate is unambiguous. The 

D 

noiseless coarse correlation result L c [0]=we v 71 co =we D , which 
simultaneously represents the expected value of the correlation, and has the phase 
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J^syncr i_. r ^ D 



angle 2rc ^ c ^ f - Clearly, 5/ < — — must be provided in the noiseless 



D 2D svnc 

sync 



case to ensure non-ambiguity of the coarse estimate. If noise is present, a 
synchronization failure rate will occur, depending on the difference between the 
actual £y-and this upper limit. 

This coarse ambiguity problem occurs for large ^only, but apart 
from that the coarse frequency offset estimate obtained from L c [0] must 
additionally lie within the unambiguity range of the fine estimator. This will turn 
out to be the more stringent problem. Thus the coarse estimator must at least 
provide that the residual frequency offset after the first correction lies such that 
|27tA£ o fc 0 71 < k , i.e., the fine estimator may work unambiguously. This phase 

D 

limit corresponds to an argument of ±tc — which defines the decision regions 



in the coarse correlation result. Consequently, the phase range for arg(L [0]), 
where the fine frequency synchronization does not fail is given by 



arg(L%0]y2iz^Z f 



D 

< 7t g??g 



FIG. 7 is a phasor diagram of the noiseless 



(coarse) correlation L c [0] (expected value) and the non-ambiguity region of the 
fine estimator for the sandamble structure. The phase range where the fine 
frequency synchronization does not fail is depicted as non-hatched area in FIG. 7, 

and as already indicated there, the opening angle tc — ^ of the decision region 



is usually smaller than tt/2 for a reasonable choice of OFDM parameters. 
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Now, we have to distinguish two cases to evaluate the probability of 
false lock due to one of the reasons mentioned above: 

The non-hatched area in FIG. 7 lies entirely within the non- 
ambiguity range of the coarse estimator, i.e., entirely above or below the negative 

D D 

5 real axis. This is mathematically expressed by n-2n ^ wc l^| > rc— 22£ , i.e., 



|^| < — ~ — . In this case we have symmetrical conditions for the 

2D sync k 0 



error events with respect to both decision region boundaries. They are 
overstepped if the magnitude of the boundary-orthogonal noise component - with 

2 D 
variance wo n (per real dimension) - exceeds the amplitude value irsinTi ^ 

k 0 



10 (see FIG. 7), in either one of the two directions. 

The second case is valid for large relative frequency offsets, i.e., in 

the range — ~ - — — < |£J < — — — , where the distances of the, noiseless 

7D k f 2D 



(expected) correlation point L c [0] to the two failure boundaries are differing 



and given by trsinii 



f 2D ) 
V u J 



and u'smn—^ f respectively. 



15 



The probability of frequency synchronization failures, P ff , in the 
coarse/fine scheme can therefore be approximated as 
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2-P 



chi 



D , 



v o 



D K-D 



sync 



P ~ 4 



ck 



sinTt-p.oJ 



.srync U 



(35) 



sinTi 



2D 
1— f=Kyl 



Z> k Q- D sync^, D 



I 



2D 



sync 



where P ch i(tf,fc) is given by: 



IN 



— w 



(aJ) D ^ c -l)! 

D -1 



Z) -1 

w c exp 



!+_££_ 



2 2 



1 — 



1 + 



2b 



(35.1) 



and has to be evaluated for calculation. 

Note that we assumed for the presented derivation, that the 
frequency synchronization failure probability is merely caused by the failure of the 
coarse estimator. This assumption is usually satisfied, because the estimation 
variance of the fine frequency offset estimator is orders of magnitude smaller than 
that of the coarse estimator. More specifically, the frequency estimate obtained 
from the fine correlation result can be interpreted as noiseless (near-zero variance), 
when compared to the coarse estimate. 

Second Technique: Fine/Discrete (TP) Approach 
In accordance with an aspect of the invention, the frequency offset 
is first estimated with the fine estimator and then corrected (even though eventually 
with discrete-valued ambiguity). An attempt is then made to finally resolve this 
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ambiguity by a postprocessing unit operating on one or more correlation results 

with a different correlation basis. The failure rate performance in the primitive 

version of this FD approach (described more fully below with reference to Eq. 

(37)) can be comparable or even identical to that of the CF, but it offers some 

5 further advantageous features, which should be exploited to significantly improve 

its performance over that of CF. 

According to the preamble structure given in FIG. 4, the fine 

(2) 

estimator uses the correlation basis kQ—D^c+D+Dg to estimate the frequency 
offset with very low variance. The estimate exhibits a normalized ambiguity 

10 interval of D/k 0 . This means that after "correction" of the estimated fine 

frequency offset, the received signal can be interpreted to exhibit solely uniformly- 
spaced and discrete-valued frequency offsets. The operating principle of the fine 
correlation and correction unit is illustrated in FIG. 8. In the following notation, 
A/ co and £y denote the remaining offset values after the already described 

15 frequency offset correction according to the sandamble has been accomplished. 

Thus the remaining discrete- valued frequency offset is determined by the property 
2nAf co Tk 0 = 2nA f with A^eZ. This results in 

Afco = Tf A f or ^ = f A ' ' A ^ Z < 36 > 

where Ay is the integer- valued parameter, indicating the ambiguity interval, the 
number of which must be detected by the new postprocessing stage. 

20 In FIG. 8, the horizontal axis represents the relative frequency 

offsets ^-and the axis is subdivided into the normalized ambiguity intervals of 
length D/k 0 . The fine/discrete algorithm works such that the original continuous- 
valued ^somewhere on this axis is estimated by the first (fine) correlation stage. 
Due to the ambiguity, this normalized offset estimate will be equal to the distance 

25 relative to the nearest interval centroid and therefore not necessarily equal to In 
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FIG. 8 the crosses indicate some continuous-valued f;y and it is shown to which 
centroid (bullets) they are shifted ("corrected") by the first processing stage. 

Now, if it were possible to perfectly estimate (or rather detect) the 
frequency offset correctly without limitation on the absolute magnitude of ^or at 
5 least for an acceptable number of intervals around 0, a frequency estimator with 

a variance identical to that of the fine estimate but with a large lock-in range would 
result. 

For this purpose the same correlation as in Eq. (29) is used, where 
the samples r[k] now denote the received signal samples after the first (fine) 
10 frequency correction. In the FD case, we obtain for the simple correlation 

L°[k] = L°\k] with 

L a [k] = £ r * [k - K y[k+D -k]. (37) 



15 



D c is typically in the range of about 10 samples to be correlated. This is the 
correlation performed by the second correlation unit 515 shown in FIG. 5. The 
principal position of the samples exploited in the correlation is depicted in FIG. 6. 

As in CF we assume perfect frame synchronization, and therefore 
we find in analogy to Eq. (31) and with Eq. (36). 



£ Cl [0]=E |r[-K]| 2 e +j27tA/ - D -^£ | rT -K]| 2 e + ^^A f 



(38) 



for the useful correlation component at the correct timing instant (k=0). Thus the 

cl 

expected argument of the correlation L [0] is 



arg 



t~ cl [0]) = 



' D N 
k n 1 



mod2ir, 



(39) 
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where the modulo operation here denotes the reduction to the argument (phase) 
interval [-*rc, ir). Considering the expression in Eq. (39), we find that, for suitable 
values of Dsync , the correlation result can take distinguishable phase angles for 

K 

values, A f , of practical interest. FIG. 9 depicts a Table 900 that presents an 
5 overview of the phase angle values associated with the frequency offset centroid in 

the interval A f for various OFDM parameter values. More specifically, the Table 

900 presents a summary of discrete circle fractions — arg(£ d [0]) e [-0.5, 0.5) 

2n 

(phase angle normalized to 2n) associated with the centroid of interval A f for 

(2) 

various OFDM parameters. The fixed parameters are £> g = D g = 8. The value 
10 D/kQ gives an idea of which NCFO ^(normalized to the respective D) is 

associated with the discrete frequency offset interval A f = 1 . The latter value is 

identical with the centroid spacing. 

It is obvious from the Table 900 that there are more or less 

advantageous parameter combinations and that in the specific case the appropriate 
15 parameter ratio, — , has to be designed to meet the needs of the maximally 

occurring oscillator instabilities. From the parameter set in the first table row 

H.=— I we learn that in this case an estimation is only possible for | A f | <3, 

because the circle is uniformly subdivided into only seven detection intervals. 
This gives a maximum normalized frequency estimation range of 
D 



=2, which is quite a lot when compared to the specified 10 ppm 



oscillator instabilities. In an exemplary system in which D = 64 and/c = 5.2 
GHz, this results in max ^ ~ 0.266. Without the newly proposed discrete 
estimator stage (postprocessing), the range of the sandamble structure would be 
restricted to |^|<— •— =0.286. This would actually be just enough for the 
25 exemplary system (with the same OFDM parameters), so that with the specified 10 
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ppm the discrete stage is not required. If it is later desired to use cheaper 
oscillators, this discrete stage might become necessary. 

Another point to be mentioned in favor of FD is the fact that 
depending on the ratio -2. 9 the maximum achievable estimation (lock-in) range 
5 may even become larger than that of the CF approach. In CF it is inherently 

restricted to |£y|< — — — =1 for the above discussed parameters D = 64 and 

^^sync 

For the other parameter choices in FD, no double representation 
occurs for the whole range of A f given in the Table 900. However, the discrete 
10 points clearly become more dense on the circle. The parameter set in the second 



15 



row 



< D ^ 



provides an exactly uniform partitioning of the circle into nine 



discrete frequency points (i.e., decision regions). For | A^| ^5 » ambiguity will 
occur there as well, so that the maximum normalized frequency estimation range is 



again bounded by |^| 



— +4 



•-—-=2 . For the parameters in the lower row 



k~2\ 



, the range of permitted Af before ambiguity occurs is larger than the 



±5 listed in the Table 900. The lock-in range is actually ± 10 with very dense but 
uniformly spaced points on the circle for frequency decision, and thus we have 

| £y| <^ = 8 • But clearly, this range must necessarily be larger in realistic 

(i.e., practical) systems, because with D = 256 subcarriers and the same absolute 
20 frequency offset A/ co . (i.e., the same RF front ends) the NCFO ^-can become four 

times that large, when compared to the 64-carrier OFDM system. The FD system 
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in the lower row could theoretically cover the same range of absolute frequency 
oscillator instabilities as the parameter sets in the first two rows. 

Once again it is mentioned that the first row of the Table 900 
actually represents a uniform partitioning into 7 discrete points if all values | A f | s3 
5 were to be used in the detection process. For illustration and further analysis, this 

phasor diagram is plotted in FIG. 10 for a relaxed detection up to | A f | <2 only 
(i.e., for the five points: -2, -1, 0, 1 and 2). This still provides a lock-in range of 

|£J<(0.5+2)-— = 1.428, which is more than sufficient for the parameters of the 

exemplary system mentioned above. With such a FD frequency estimator the 

10 carrier frequency oscillator accuracy requirement could be relaxed to over 50 ppm. 

Consequently, the correlation result £ ci^ has to be used for the 

estimation of the discrete frequency offset. The performance criterion of interest 

is clearly the error rate in this discrete estimation (or, more appropriately, 

detection) problem. Because the exploited correlation is the same as in CF, the 

15 statistical properties of £ cl ^ are identical to the ones of £ c [q] • Hence, the 

derivation in CF applies here as well. 

As an example, we analyze the performance for the case of the 

(2) 

correlation constellation depicted in FIG. 10 for D g = D g v = 8, D = 64 and 

D sync = 32 and | A f | ^2. We assume that the discrete frequency offset lies in 

20 interval zero (i.e., Af = 0) and calculate the probability that the discrete frequency 

estimate is incorrect. The probability of error can be approximated by neglecting 

the points with | A f | = 2 and applying the union bound yielding two times the 

probability that the orthogonal noise component in ^ cij-Qj is larger than half the 

distance between the correlation signal points A f = 0 and A f = 1. This half 

25 distance is given by u ■ sin it (see FIG. 10). Similar calculations have to be 

K 

performed to obtain false-lock rates for the other discrete frequency offset intervals 
Af* 0. 
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The probability P ff of frequency synchronization failure in the 
specific constellation diagram of FIG. 10 can therefore be approximated as 



K o 



D 



ff I chl Ir n chlV 2k n 'lb f IJr f 



^0 



2k 0 * 



3 — < 1 5 J <— or | A f | =2 



(40) 



where P chi according to Eq. (35.1) has to be used for calculation. 

It is obvious that the case A f = 0 results in the lowest error rates, 
5 while false-lock probabilities for higher offsets are larger. This estimator property 

is quite reasonable. 

As a first conclusion of all this theory, we can compare the false- 
lock performance of the CF estimator according to Eq. (35) and the FD approach 
from Eq. (40) and we find that for lower frequency offsets (i.e., \ty<-^- ), the 
10 error rate is identical in both schemes for this special choice of OFDM parameters. 

Third Technique: Improving the Power Efficiency of the Detection 
in FD 

So far, it has been shown analytically that for the parameter choice 
D = 64 and D sync = 32 in the sandamble scheme, the FD approach with | A f | <; 2 
15 is equivalent to the CF approach in terms of false-lock rate at low frequency offset 

values. 

To mention a first advantage of the FD approach, it is pointed out 
that a second correlation can be exploited to determine the discrete frequency 
offset, so that a slight gain in power efficiency of the detector can be achieved. 
20 For this purpose, the second correlation 
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L c2 [k] = E^[^K]r[^ + Z) sync ^ 0 -K] (41) 



is examined more closely. This is the correlation between the G 2 and B 2 regions 

that is performed by the third correlation unit 517 (see FIG. 5a). The principal 

c2 

sample positions to be used to calculate L [k] have already been indicated in FIG. 
6. For the optimum time instant (k = 0), it follows in analogy with Eq. (31) and 
5 with Eq. (36) that the noiseless component reads 

L C \0] = £ kT-*]| 2 * v ' 2nAf - {D ^ T = £ \f[-K]\ 2 e k * . (42) 

*=1 K = l 



C2 r 



Hence the expected argument of L [0] is 



arg 



t~ c2 [0]) - 



( 



2tc sync ° A, 
it. 7 



mod27t = 



mod 2n, 



(43) 



which turns out to be the same expected phase as that of L [0]. Note that the 
discovered equivalence between the arguments is established only by the operation 
of the fine frequency correction stage, which performed earlier on the received 
10 samples. 

Thus the parameter to be exploited in the proposed FD approach 
with increased power efficiency is the correlation sum 

L c [k] = L cI [k] + L c2 [k] . (44) 
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In the exemplary embodiment, this correlation sum is generated by the combiner 
519 (see FIG. 5). Because optimum time synchronization has been assumed, we 
evaluate it at the perfect time position (k=0) and obtain: 

L cI [0]+L c2 [0]*J22\n-K]\ 2 e 

K = l 

+ E r [-K](«[^ nc -K] ^P^+Ko-K]) 

K = l 

E< r[D -k] 

k=i r*[-K] 



= £2|r[-K]| 2 e "> 



K=l 



sync -> sywc 0 J ^j-_ K j 



so that the useful component now is 2we 0 , with the same definition of u 
5 as previously introduced in Eq. (42). Thus, the amplitude of the useful signal 

component is doubled. After a short calculation, the variance of the noise 
component is found to be u-6o 2 n . The probability P ff of frequency 
synchronization failure in the specific constellation diagram of FIG. 10 can 
therefore be approximated as 
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2-P ch ^m^,3o 2 n ), \i f \ A f =0 

P^2smn^3o^ + P^2sin^,3olx ^<|^<|^or |A f | =1 , (46) 
2^(23^^,30^ || <| 5/ | <|| or IM =2 



where again P chi (a, &) according to Eq. (35.1) has to be used for calculation. 

From the comparison of the failure rates analytically obtained from 
Eq. (40) and Eq. (46), one can conclude that for any £y the same false-lock 
performance is now achieved with a receive power reduced by approximately 
1.25 dB, when compared to the single correlator ("primitive" version). 
Consequently, the use of the additional correlation slightly improves the power 
efficiency of the discrete frequency estimator. 

FIG. 1 1 is a graph in which the theoretical failure rate performance 
of CF is compared to that of FD over a range of E s /N 0 values and over a range of 
relative offset Rvalues. For both approaches, the guard interval of the first 

synchronization symbol is extended to = 14 , so that a (quite low) number 

of D c = 6 correlation products is used in the coarse correlation which can be 
interpreted as a diversity factor of 6 for the detection problem. 

FIG. 12 is a graph in which the same failure rate comparison is 

15 provided for = 18 , which implies, that D c = 10 correlation products are 

used in the coarse correlation. It should become clear that, due to the higher 
diversity factor, the performance is improved. The curves for the false-lock rates 
now decrease much steeper for rising SNR. 



5 



10 



WO 99/65180 



PCT/EP99/03793 



-56- 

The diagrams show that the CF approach exhibits a constant 
performance over a wide range of ^up to \J£ f \ < — — — ^ =0.714. 

Increasing the frequency offset beyond this limit results in a considerable and 
continuous performance degradation. The complete break-down of CF is reached 

5 at | £y| =1 • But clearly the noise performance of CF is already unacceptable 

beyond =0.9 , so that the full nominal lock-in range of \^ f \<\ is only 

valid for noiseless systems. 

The performance characteristic of the FD approach is markedly 
different, as the performance is constant within one interval, but degrades (or 

10 improves) in hard steps. If it were necessary to estimate the maximum possible 

range | A f | <3, the error rate in all intervals would be identical to that for | A f | = 2 
(diamonds for FD in FIGS. 11 and 12), because then the circle is uniformly 
partitioned in this case. It is worth noting that if it were only necessary to estimate 
| Ay | < 1, the false-lock rate for | A f | =1 would be even lower than the one for 

15 A f = 0, so that FD would outperform CF in the entire range up to | £ f | = 0.857. 

Another point which should further be considered is the following: 
So far we have only considered preamble structures in which 

£> (2) ~ D and D^ l) >D , that is, those preamble structures in which the first 

part of the split synchronization preamble provided the additional samples for the 
20 coarse correlation. Clearly, this is not compulsory because, in alternative 

embodiments, the second part of the synchronization preamble could equivalently 
be transmitted with the extended guard interval, while the first part could be 
transmitted using the conventional guard interval duration. This would affect the 
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results so-far presented in this section in several ways. First, the correlation basis 
ko of the fine correlation would become adjustable in certain limits, but apart from 
this parameter nothing would be changed in the algorithm during the first (fine) 
estimation stage. It is worthy of note that the ratio 

— = — which determines the FD decision regions for the 

u sync g g 

discrete frequency detection unit would be adjustable (in limits) through the 
parameter k 0 or rather by D^ 2) . The most important thing to point out is that the 

correlation sum for the coarse estimate would differ considerably as we now 
exploit 

L cl [k] = E^[* + *b-KW*+*b+^-K] (47) 

K=l 



10 and 

L c2 [k] = £r*[*+Z> -K]r[*+V K ]- (48) 
The temporal position of the associated correlation windows which are used to 

1 r> r ) 

calculate L [k] and L [k] when this alternative sandamble structure 1301 is used 
are illustrated in FIG. 13. It can be seen that an extended guard interval occurs in 
the second part. 

1 5 For the two single noiseless correlation components - which 

simultaneously represent the expected correlation values - we find now at the 
correct timing instant 
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K = l 



(49) 



K = l 



L C \0] = E^[^"K]r [*b-K]=E l^[>t 0 -K]| 2 e k « . (50) 

K = l K = l 



Therefore, the expected argument values of the two correlation results in this case 
(D?>D?) are 



arg (r'[0]) = 



< D X 

k n 1 



mod 2ic 



(51) 



arg 



t~ c2 to]) = 



2-k- 



K-D. 



A* 7 



mod 2tc = - 



Jfc. 7 



mod2n;. 



(52) 



For the first correlation we obtain the same mean argument as 
before, while the second correlation now exhibits the exactly negated argument. 
Thus the correlation sum to be exploited to achieve maximum detection power 

efficiency in the case D g (2) >D^ l) reads 

L c [k] = L cI [k]+{L c2 [k])\ (53) 
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which has the same statistical features as the "double correlation" in Eq. (44) for 
the case Dg l) >Dg 2) . Therefore, both variants achieve nearly the same 

performance if the slightly increased fine correlation basis (and therefore 
marginally lower estimation variance and negligibly reduced lock-in range) in the 
5 latter approach is neglected. It is noted that in this embodiment, the simple 

combiner 519 depicted in FIG. 5a would have to be supplemented with a complex 
conjugation unit (not shown) at one of its inputs, in order to receive the value 

L c2 [k] ' 

Usin g A Priori Probabilities for the Detection in FD 
10 So far the discussion has not considered the pdf of the frequency 

offsets which will actually occur during practical operation, and clearly they will 
not be uniformly distributed. It is justifiable to assume that NCFOs around zero 
will be predominant, while offsets in intervals farther out will be rather seldom. 
In the simplest model of uniformly distributed and statistically independent 
15 frequency offsets in the transmitter and the receiver oscillator, the overall 

frequency offset pdf will follow a triangular shape. Thus we are not interested in 
minimizing the probability of frequency failure for each of the intervals, but we 
would like to minimize the overall a posteriori frequency failure rate by taking into 
account the pdf of NCFOs, that is, considering the a priori probabilities of 
20 estimation intervals. The decision rule to be applied can be interpreted as that of a 

digital transmission with a generally non-uniform PSK signal set. The correlation 
phase angles (signal points) are defined by the respective correlation results and - 
most importantly - occur with unequal a priori probabilities. 

The decision A f on the discrete frequency interval shall be in 

25 favor of the a posteriori most probable one. This is expressed in the decision rule 
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A f = argmax Pr {A^ = A f \L c | , 



(54) 



where for brevity of notation we introduced L c = L c [0] for the correlation result 

at the optimum timing instant. Note that here L°[k] stands either for L° l [k], L° 2 [k] 
or the suitable combination of both. 

Now we focus on the most power efficient case of the double 

correlation, where the useful component was 2u -e ^ and the Gaussian 



2 

noise variance was u - 6o n (per complex dimension). We modify the a 

posteriori probability in several steps and obtain the expression 
Pr{A f = A f |Z c } 

P L <( L ) 

- vJhlA . _L_ . exp ( -, Z <-2„e ^\/[6u ■ A (55) 



gjgf - A f} 
6nuo n - p L c(L c ) 



Dsync - A 



exp^-|Z c | 2 + 4w3t{e k ° ' ■ i c } -4u 2 j /(6u • 



where Pr{A f =A f } is the a priori frequency offset interval probability which 

can be obtained by integration of the frequency offset pdf within the respective 
10 interval limits. The interval -conditioned pdf P L c(L c \ A f =A f ) is in very good 

approximation Gaussian at higher SNRs, where the noise-noise products are 
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negligible (see Eq. (31)). Neglecting the constant-valued multiplicative factors, 
which do not affect the maximum probability decision because they are common 

for all we yield the simplified probabilistic decision rule 



A f = argmax 



Pr{A f = A f } - exp [l8t{e ** * • L*}J? - oj)) 



(56) 



which is independent of the centrally chi-square distributed random variable u. It 
is this decision rule that is preferably applied by the frequency offset interval 
detection unit 521 (see FIG. 5a). Equivalently, the decision rule in the logarithmic 
domain reads 



A f = argmax 

A, 



-9t{e 



-~ Dsync r 
-j27t— j- — A f 



L c ) + lnPr{A f - A f | 



(57) 



Thus the knowledge of the noise variance o n as well as the a priori 

probabilities Pr{A f = A f } of the discrete frequency offset intervals A f are 

10 required for an optimum decision on the discrete frequency offset. 

For a priori equally probable offset intervals, the decision rule 
simplifies to the well known type of 

A f = argmax 9t{e *° • L c }. (58) 
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Without explicit derivation we give the final decision rule for the less power 
efficient single correlation case, which reads 



A f = argmax 
if 



1 -P^if 

— 3t{e *° - L c } + lnPr{A f = A f } 



(59) 



reflecting the loss of io log 1—1 = - 1 25 dB m power efficiency between single 

6iol 4 l • 

and double correlation. This follows directly from the factor -± in the amplitude 

3 

factor, when comparing Eq. (57) to Eq. (59) and the fact that the correlation result 
L° in Eq. (57) will be twice as large for the double correlation case. 



Simulation Results with Perfect Time Synchronization 
First, simulation results for the CF scheme in AWGN (non- 
dispersive channel) and a multipath channel according to a decaying delay-power 

10 profile 1 are presented and compared to the theoretical results obtained from the 

AWGN performance analysis. The simulation parameters for all presented results 
utilize D — 64 and D u = 53. 

FIG. 14 is a graph in which the false-lock rate P ff of CF in an 
AWGN channel is shown for various normalized frequency offsets |£ f | e [0.0,0.9]. 

15 Asymptotically, the performance is exactly predicted by Eq. (35). Even the 

predicted continuous degradation for |£ f | k 0.714 is confirmed by the simulation 
results. For lower SNRs, there is a slight difference between the simulation results 
and the analytic results which is due to the neglected noise-noise products in the 
theoretical derivation following the steps from (29) to (31). 

20 FIG. 15 is a graph in which the frequency synchronization 

performance of CF for transmission over a severe multipath channel is shown. 
The theoretical predictions incorporated in the plot are valid for AWGN and 
surprisingly, the performance difference is very low: approximately 1 dB. The 
difference is caused by the channel-induced correlation of subsequently received 
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samples. Hence the D c samples in the correlation sum are no longer statistically 
independent, so that the useful correlation component does not follow the chi- 
square distribution which has been assumed in the derivation. 

Thus far, we can conclude that all assumptions in the analytical 
5 derivation of the false-lock rates for AWGN were justified and the performance 

results are verified via the simulated CF synchronization behavior. Even though 
the analysis was performed for AWGN, the analytic result can be used to predict 
the performance for severe multipath channels within an uncertainty of roughly 
1 dB or less. 

10 What follows are simulation results for the FD scheme with "double 

correlation", able to detect 5 discrete intervals, i.e., | A f | <; 2. This gives a lock-in 
range of | £ f | <; 1.429, which is 43 % larger than the nominal range of CF. It 
should be noted that CF actually cannot exploit the total nominal lock-in range, 
because the performance is continuously degrading up to the total break-down at 

15 the nominal limit | £ f | = 1.0. FIG. 16 is a graph that depicts the false 

synchronization probability for transmission over the AWGN channel for various 
values of £ f . The step-wise degradation of the frequency synchronization behavior 
is entirely confirmed by the simulation results. Asymptotically, the theoretical 
false-lock rate is nearly identical with the simulated results, and the difference at 

20 lower SNRs has its origin (like in CF) in the neglected noise-noise products in the 

correlation sum (see Eq. (44)). 

In a multipath channel, the FD approach shows a good 
correspondence to the theoretical prediction derived under the assumption of a 
transmission over AWGN. This performance evaluation is shown in the graph of 

25 FIG. 17. As in the CF scheme, the multipath simulation results are within 1 dB 

from the AWGN prediction, so that the theoretical result from Eq. (46) serves as a 
good mean to evaluate the frequency synchronization failure rates even in severe 
multipath channels. 
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To conclude these simulation results, the benefit of using a priori 
information for the occurring frequency offsets in the FD detection rule is 
motivated via the simulation result in the graph of FIG. 18, which has been 
obtained on an AWGN channel. The parameter setting was such that the a priori 
probability for the zeroth frequency offset interval (A f = 0) was 0.90, while the 
adjacent two intervals | A f | = 1 and | A f | = 2 obtained a probability of 0.04 and 
0.01, respectively. It is apparent that the false-lock rate P ff is drastically reduced 
for interval zero and slightly reduced for | A f | = 1. The costs for this reduction are 
to be seen in the increase of P ff for the intervals with | A f | = 2. Hence, the 
communication links with large frequency offsets and quite low probability will 
suffer a higher false sync rate. The overall synchronization failure rate is 
minimized by this detection rule. 

Fine/Discrete Approach with Pilot Corre lation (FDPO 
An approach is now described that is not restricted to application 
only in the FD approach with the sandamble structure as described earlier. 
Rather, this frequency-domain correlation can as well be used in any approach 
with sufficiently high discrete-valued frequency ambiguity, i.e., |£ f | ^ 0.5 or even 
integer values of the residual NCFO £ f . It is only required that a frequency 
ambiguity with an interval known at the receiver be generated by some fine 
frequency correction stage. 

The core of this Fine/Discrete approach with Pilot Correlation 
(FDPC) is the exploitation of differential pilot symbols across the subcarriers of 
the frequency multiplex. These may be the same pilots as will be described for use 
in an (iterative) frequency-domain frame synchronization aspect of the invention. 
Thus this approach does not require any additional redundancy. The received and 
fine frequency corrected signal is frequency shifted (mostly in the time-domain by 
modulation) with integer multiples of the frequency ambiguity interval and 
transformed into the frequency domain by a separate DFT each. If the relative 
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frequency offset associated with the frequency ambiguity interval is an integer 
multiple of a subcarrier spacing, this frequency shift can be done directly in the 
frequency domain by a shift of the subcarrier amplitudes. Hence only one DFT 
would be required in this case. 
5 After transformation into the frequency domain samples y the 

correlation, L, is determined for each of these trial signals, / v , in accordance 
with: 

L d - Yi r M 7 M+ £v * (60) 

Vve/ 



It will be seen that this same frequency-domain correlation is useful in the stage 2 
aspect of the invention, to be described later in this disclosure. 

10 Now, instead of evaluating the phase angle, the amplitude of the 

correlation result L is evaluated and the interval which yields the maximum 
correlation amplitude is selected. 

It is now advantageous to restrict the single observation blocks in 
the estimator to K = 2, k$ = 1, (2kQ = K) i.e., D cig = 1 to keep the radial noise 

15 component low. Unfortunately, this choice has been found to be suboptimum for 

the frame synchronization algorithm, because there the phase is evaluated. For 
phase estimation, it is desirable to use D cig > 1 . Hence a compromise in pilot 
arrangement between the frame synchronizer and the discrete frequency detector 
needs to be found, based on the specific problem setting. 

20 Apart from that, this approach is data limited, and will therefore 

suffer a flattening, even in noise- and interference-free situations. This means that 
owing to the random data within the OFDM symbol which carries the pilots, a 
false lock may be caused by eventual match of random data with the correlation 
pattern. The error floor can be moved to lower values by increasing D cg , so that 

25 the probability of random data match gets very low. If an OFDM sync preamble 

without random data is used (i.e., it actually is a training symbol in the strict 
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sense), the subcarrier values in this symbol can be optimized to prevent the 
flattening. 

Simulation Results with Imperfect Time Synchronization 
FIG. 19 is a graph depicting the frequency false-lock rate P ff of the 
5 CF scheme with D c = 6 and various misplacements D m in a severe multipath 

channel. Various normalized initial frequency offsets |£ f | e [0.0,0.8] are used for 
simulation. The time sync error causes an error floor for P ff . 

The performance of the FD approach is given in the graph of FIG. 
20 for equal a priori probabilities of the frequency offset intervals, while the graph 
10 depicted in FIG. 21 reflects the performance characteristic with the same a priori 

parameters previously used. It can be seen from a comparison of FIG. 20 with 
FIG. 19 that at low frequency offsets FD is less affected by frame synchronization 
errors than CF. Including FIG. 21 in the comparison, it can be seen that this 
robustness becomes increasingly large for unequal frequency offset intervals. It 
15 should be noted that the length of the prolonged guard interval of D c = 6 is quite 

low, when compared to the misplacements D m . An increase in D c would lower 
the error floors. 

Simulation results for the FDPC approach are now presented. From 
the graph depicted in FIG. 22, it can be seen that this approach produces 
20 unsatisfactory results in frequency -selective channels. Only D cg = 10 together 

with D r :„ = 1 were used in this simulation, which might be too few pilots for a 

Clg 

. reliable decision in a multipath channel. Considering the graph depicted in FIG. 
23, we learn that the performance of FDPC in AWGN is not that bad for an 
approach that exploits pre-existing redundancy, that is, an approach that does not 
25 introduce D c additional training data samples (prolonged guard interval) as CF or 

FD requires. A further increase of D cg might yield better results in AWGN as 
well as in multipath channels. 
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Alternative Embodiment of First Stage 

An alternative embodiment of the first stage will now be described 
with respect to FIG. 5b. From a mathematical point of view, this embodiment 
accomplishes the same thing as the embodiment described above with respect to 
5 FIG. 5a. However, this alternative embodiment is less expensive because it 

requires fewer multiplications and less memory. In this embodiment, received 
(noisy) samples r[k] are supplied to a sample memory 501. It is known that the 
stored received samples r[k] will include a sandamble structure 503 . In a first step 
of the first stage, the received samples r[k] are supplied to a first correlation unit 

10 505 that generates a coarse frame start (time) estimate 509 as well as a first 

frequency correction 507. These are generated by determining similarity metrics 
from the received samples r[£], and using information obtained from the generated 
metrics to detect the sandamble 503. From the detected sandamble 503, both the 
first frequency correction 507 and the coarse frame start (time) estimate 509 may 

15 be determined. Alternative embodiments of the invention may be derived by using 

different types of similarity metrics, as explained above with respect to FIG. 5a. 

In order to obtain a decision variable, G 2 is correlated with A 2 by a 
second correlation unit 515. In addition, a third correlation unit 517 correlates the 
G 2 region with the B 2 region. In order for the two correlation results to be 

20 combined, it is necessary to perform a frequency correction on each. In this 

embodiment, this is accomplished by respective first frequency correction units 
518. Each of the first frequency correction units 518 receives a correlation result 
from a respective one of the second and third correlation units 515, 518, and 
corrects this as a function of the first frequency correction 507 generated earlier by 

25 the first correlation unit 505. The frequency correction units 518 operate by 

multiplying the received correlation value by a complex-valued phase rotating 
factor which is calculated from the first frequency correction 507. After each has 
been frequency corrected, the two correlation results (from the second and third 
correlation units 515, 517) are combined (e.g., by addition) by a combiner 519, 
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although such combination is not essential to the invention. That is, in alternative 
embodiments, only one or the other of the two mentioned correlations need be 
performed. Returning now to the exemplary embodiment in which both 
correlation results are generated and combined, the combined correlation result is 
5 supplied to a frequency offset interval detection unit 521, which detects a second 

frequency offset 523. This frequency interval decision 523 is combined with the 
first frequency correction estimate 507 from the first correlation unit 505 to form 
an overall frequency offset estimate 525 (e.g., by multiplying the stored sample by 
a complex- valued rotation factor sequence that is a function of the estimated 
10 frequency offset value). The overall frequency offset estimate 525 may then be 

applied to frequency correct the received samples r[k] . Other aspects of the 
particular operations performed by this alternative embodiment are the same as 
those described above with reference to FIG. 5a. 

SECOND STAGE: FINE FRAME SYNCHRONIZATION FOR OFDM 

In this section a post-processing stage for frame synchronization is 
outlined that achieves a more accurate estimate in packetized spontaneous OFDM 
transmission. To estimate the time offset of the received modulation signal a ML 
timing offset estimation from the frequency-domain is chosen that applies a 
correlation technique on frequency-domain received subcarrier amplitudes. The 
algorithm starts operation upon a coarse timing estimate obtained from a time- 
domain metric. Other techniques might detect Null symbols, such as is described 
in Markku Kiviranta and Aarne Mammela, "Coarse Frame Synchronisation 
Structures in OFDM", ACTS Mobile Telecommunications Summit , pages 464- 
470, Granada Spain, 1996 (henceforth, "[KM96]"), or simply the maximum of a 
power detector. 

The performance characteristic of this frequency-domain post 
processing unit is derived in the following for an application to the synchronization 



15 



20 



25 
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symbol with £> sync carriers, but it should be kept in mind, that this post processing 
is applicable on each regular Z)-carrier symbol as well. 

The Frequency-Domain Frame Synchronization Approach 
We assume perfect frequency synchronization and find 

Y M =A MV H[v]e aync < 61 > 



5 for the noiseless receive subcarrier amplitudes, f > in subcarrier v if the 

At 

Normalized (sampling) Time Offset (NTO), £ r = , is not too large, a is 

T ** ,v 

the transmitted subcarrier value in subcarrier v, and H[v] ls the channel influence. 
Hence a time offset has basically the effect of a linear phase term across the 
received subcarrier amplitudes. The benefit of exploiting the fundamental system 

10 theoretical property in Eq. (61) for time synchronization is mentioned in S.U. 

Zaman and K.W. Yates, "Multitone Synchronization for Fading Channels", 
Proceedings of the International Conference on Communications (ICC'94) . pages 
946-949, New Orleans, USA, 1994 (henceforth, "[ZY94]*). There, the authors 
try to apply linear regression methods to find axis crossings and slope of the linear 

15 increasing phase rotation across the subcarrier multiplex. From these quantities 

they derive estimates for the carrier phase and time offset. The solution now 
described in this disclosure works differently. We interpret (some of) the 
subcarrier transitions from subcarrier amplitude v to ^ >v +l t0 ^ e differentially 
encoded in the frequency direction with the complex-valued information amplitude 

20 / M v (differential pilot) so that 

A M.v.r A M.v mI » v Vv e 1 (62) 
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with the subcarrier index set 

/«{vO^Vl'^.^-l;y=0 > ...^ cfc -l} c [O.D-2]. (63) 

This set expresses thatZ? cg ("eg" means "correlation groups") 
blocks with D cig ("cig" means "correlations in group") directly adjacent 
differential information symbols / v (are interpreted to) perform a differential 
5 encoding across the subcarriers. The single correlation blocks are assumed to be 

spaced D cs > D cig + 1 subcarriers apart. The principle arrangement of these 
differential subcarrier pilots within the subcarrier multiplex is depicted in FIG. 24. 

Considering Eq. (61), it is apparent that the correlation of 
frequentially adjacent subcarrier amplitudes can provide some reasonable means 
10 for estimating the relative time offset For the noiseless received subcarrier 

amplitudes, the single correlation product equates to 

K v K =K v\ v + i H*[ v] H [ in- 1] e + D ~ 1 , < 64 > 



(where ^ * ^ ji/ ) and we observe that the linear increasing phase is 

transformed to a constant phase offset. Especially in an AWGN channel, where 
the discrete channel transfer function is H [v] = 1 , V v, the mean value of this 
15 estimate is a perfect estimate of E, r 

The dependence on the argument of the differential information 
symbol 7 AV must be canceled by multiplication with the complex conjugate of a 
(possibly perfect) estimate /^ v of this complex-valued amplitude (modulation 
removal). 
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For modulation removal, at least two well-known techniques can be 
applied. These are: a Data Aided (DA) technique, and a Decision Directed (DD) 
technique. 

In the DA technique, the receiver a priori knows the frequency- 
5 domain symbols (subcarrier amplitudes) dedicated for fine time synchronization 

purposes. Clearly, the known symbols represent additional redundancy if they are 
not used for further estimation procedures. They cannot be used for data 
transmission and therefore the number of these pilot symbols should not be 
excessively high. This approach is very robust and the only one that works if the 

10 number of carriers is too small. Apart from that, this approach offers the 

possibility to estimate frequency offsets larger than one subcarrier spacing, 
provided that the "fractional frequency offset" (i.e., the frequency offset 
component with magnitude smaller than half a subcarrier spacing) is corrected 
prior to this stage. For this purpose, the values of the received frequency-domain 

1 5 symbol sequence and some shifted versions of it are compared to the known 

symbols. The frequency shift with the best match is the estimate for the "integer 
frequency offset." 

In the Decision Directed (DD) technique, the estimator operates on 
detected modulation symbols, i.e., on decisions from the demodulator to remove 

20 the modulation. The differential encoding with / v can be virtual in that it does 

not necessarily represent information. Hence, no redundancy is needed for this 
type of estimation and, thus, the number of symbols exploited in the estimation 
process can advantageously be as large as the number of used subcarriers in the 
OFDM system. This kind of estimator gives satisfying performance only for large 

25 OFDM symbols (> 256 carriers) and moderate signal constellations (e.g., QPSK). 

Theoretically, the estimation variance of DD is usually better than that of DA 
because more symbols are exploited in the estimation process and adjacent 
subcarrier transitions are used, which further decreases the estimation variance. 
But the estimate of DD suffers not only from Gaussian noise and interference, but 
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also from the fact that wrong demodulator decisions produce feedback noise. 
Clearly, this approach is not as robust as DA because, depending on the cardinality 
of the modulation alphabet and the accuracy of the initial time offset, the 
demodulator symbol errors will be higher or lower due to the dominating 
5 interference noise. 

With modulation removal, the argument of the frequency-domain 
correlation of the actually received noisy amplitudes is 

def -r-^ + * * 

L = 2^ Y \x,v ^, v + l ^i,v (65) 
Vve/ 

where Y„ v are frequency-domain samples (obtained from a Discrete Fourier 
Transform (DFT) of noisy signal samples), and / * is a known or estimated 

10 discrete signal point (which is approximately equivalent to modulation removal). 

Assuming correct (error-free) modulation removal, and together with the 
restriction \ I jUtV \ = 1> ^ vand Q 2 =e j |^ | 2 j, and assuming that the correlation 
result L comprises a representative set of the available frequency bins so that the 
expected value of L can be calculated with C h [Vq] , which is actually determined by 

15 averaging over all frequency bins, we yield for the expected value of the 

correlation 

e{L}=D cg -D cig -o 2 A -c^ D ^'-C H [l], < 66 > 



with the frequential autocorrelation C H [ v 0 ] being given by 

def 1 D-\ 

CM = /r[v]//[(v+v 0 )modD] (66.1) 

D v =0 

where D is the number of carriers and ^[ v j is the channel transfer factor at 
subcarrier v . Note that we have to substitute D = D s there. Hence 
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£,=-i^arg(Z) = 5 l +-f=«g(C fl [lD < 67 > 



can be used to provide a reasonable estimate of the time offset. With the estimator 
in Eq. (67), a Active relative time offset of ^sync arg (C#[l]) is added owing to the 

271 

center of power delay gravity. But clearly, this delay is constant for the specific 
static channel parameters effective during one burst. 
5 It can be shown that the frame synchronization estimate found by 

this estimator is MMSE-optimum for differential demodulation in the frequency 
direction. 

The number of samples in the cyclic prefix guard interval and the 
number of samples in the cyclic postfix guard interval may be optimally adapted to 
10 an expected channel power delay profile. This is basically a numerical 

optimization problem that strongly depends on the multipath (echo) characteristic 
of the physical transmission link (e.g., radio link). A good compromise is to use 
as many samples in the cyclic posfix guard interval as needed to represent the time 
associated with the center of power delay gravity of the channel. 

15 A Lower Bound for the Frame Synchronization Variance 

The following analysis aims to derive a lower bound for the time 
offset estimation variance. As already outlined, the time offset estimator performs 

^J^arg^JW (68) 
2tt 2% 
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where (p is the phase increment between two subcarriers. Here, we assume that 
the two sync symbol parts of length D sync are processed with a 2D sync -point DFT. 
For this reason, the variance of the signal set in the synchronization symbol 
deserves some attention due to the symbol repetition (time diversity). Applying 
5 the theorem of Parseval, we yield the equation £> u>sync . Q 2 A = (2D sync ) • o*> 

owing to the fact that two time-domain signal parts of length D sync are exploited to 
gather the useful energy for D nsync ac ti ye ("used", hence the subscript "u") 
subcarriers. Therefore the following effective signal set variance (fictive) is true 
for the synchronization symbol: 

°5 ■ (69) 

u,sync 



10 For the noise components Parseval 1 s theorem reads 2D sync . a ^ = 2£> S yn C . Q 2 n » 

and thus we have the identity n 2 _ a 2 . 

N n 

The variance 0 * = e{( ^ -ej^}) 2 } °f normalized time 
synchronization error can be lower bounded via 



D 

2 -^sync 
O e = — 



^sync ^\i,sync 1 



2D. 



D n D * D 



Fallot 



4it 2 



(70) 



u,sync 



but it is worth noting that this variance is a purely theoretical lower limit. We 
1 5 assumed a perfectly aligned demodulation window without any interference power 

from preceding or subsequent OFDM symbols. This is clearly not true for the 
cases of typical operation. 

Instead of estimating the time offset from the frequency-domain of a 
repeated synchronization symbol, we can estimate the time offset from scattered 
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pilots being multiplexed into any one regular OFDM symbol of dimension D with 
D u non-zero subcarriers. In this case, we have no "time diversity" and the final 
result for the lower bound reads 

0 2 _ J_ D'D U 1 

" 4tt 2 D . D 2 E s /N 0 < 71 > 

eg ^Clg s u 



For D = 2D sync , this variance is 8 times larger than the variance 
5 achieved with Eq. (4.62). In the case that D = D sync we still have a 2 times 

higher estimation variance which directly reflects the diversity factor. 

Iterative Frequency-Domain OFDM Frame Synchronization 
The task to be fulfilled by the frame or time synchronization unit is 
to reliably find the starting position of an OFDM symbol in a sequence of noisy 

10 received samples that have been corrupted by ISI and additive noise on the radio 

link. As used here, the word "reliable" means that the variance of the timing 
estimate should be minimum or that the estimate is within a tolerable limit of some 
reasonable estimation mean. Furthermore, the false-lock probability should be 
minimized. In contrast to single-carrier modulation, the requirements to be met by 

15 the frame synchronizer are not that stringent, as offsets up to a few samples can be 

tolerated in OFDM. This is not the case for single-carrier schemes, where frame 
synchronization is a hit-or-miss problem. Throughout this section and for the 
simulation results that follow, we assume the worst case situation in which the 
sync symbol is embedded in a stream of OFDM symbols and is therefore not 

20 surrounded by mere noise. 
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Motivation of the Iterative Approach 

A very important statement is concerned with the interference terms 
due to the time offset in the initial (coarse) time offset estimate. This interference 
from preceding and succeeding OFDM symbols represents noise with the same 
5 power as the useful signal. It is acquired by the misalignment of the demodulation 

(DFT) window and clearly limits the estimation accuracy in the first run of the 
proposed time synchronization refinement procedure. Consequently, in 
accordance with another aspect of the invention, the estimation described above 
should be performed iteratively. The first run will usually provide an improved 

10 time offset estimate, so that a repeated realignment of the time window and a 

further run of the refinement procedure will yield increasingly improved time 
synchronization estimates with decreasing estimation variance. After a sufficient 
number of iterations and under the restriction of an initial estimate which ensures 
convergence of the iteration process, the variance of the time offset estimator 

15 approaches the theoretical lower limit in Eq. (70). Simulations for DA indicate 

that two or three runs are sufficient to approach the theoretical lower limit for 
moderate initial time offsets. 

Implementation Example 

For operation over a multipath (LSI) channel, the guard interval is 
20 advantageously split into prefix and postfix parts because the proposed estimator 

defined by Eq. (67) will add a time associated with the center of power delay 
gravity of the respective channel. More specifically, the estimated frame 
synchronization positions are rather shifted towards more positive time positions. 
Thus with the conventional OFDM systems, which have a guard interval only 
25 preceding the useful symbol part, the proposed estimator will systematically 

position the time window too far into the direction of positive time values, so that 
the demodulation window always includes interference attributable to a part of the 
guard interval of the following symbol. 
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In accordance with an aspect of the invention, this problem is 
addressed by providing an OFDM system in which a postfix is alternatively used 
alone, or together with a prefix. The format of the OFDM symbol in accordance 
with the modified scheme is illustrated in FIG. 24. As usual in OFDM, the D 
5 samples in the middle are obtained from the Inverse Discrete Fourier Transform 

(IDFT) of the subcarrier amplitudes. In accordance with an exemplary 
embodiment of the invention, a prefix D g pre is generated by a repetition of an 
ending partial sample sequence D end from the entirety, D, of the originally 
obtained samples. Similarly, a postfix D g post is generated by a repetition of a 
10 beginning partial sample sequence D begin from these D originally obtained 

samples. 

Another advantage of this structure is that the iteration will 
converge faster because, when positioned at the correct time position, there are 
guard samples both to the left and to the right of the demodulation frame, whereas 

15 in the conventional prefix structure, no guard samples will be found to the right of 

the demodulation frame. 

A block diagram of an exemplary embodiment of the stage 2 frame 
synchronization mechanism is shown in FIG. 26. This figure depicts not only the 
hardware components used in the exemplary embodiment, but also the flow of 

20 signals and order of operations, and as such also serves as a flowchart of the 

various steps performed in the embodiment. In order to depict this flow of 
operation, the same hardware component, such as the sample memory 501, may in 
some cases be depicted in several places in the figure. Notwithstanding this 
multiple depiction, only one of any such unit is required. The functions of the 

25 various units described in this and subsequent figures may be implemented in any 

of a number of ways. For example, programmable processing logic may perform 
the herein-described functions by executing program instructions that have been 
stored in a suitable computer-readable storage medium (e.g., Random Access 
Memory (RAM), magnetic storage medium, optical storage medium, and the like). 
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Alternatively, hardwired logic may be designed and built to perform some or all of 
the various functions. Each of these alternatives, as well as their equivalents, are 
intended to be considered within the scope of the invention. It will further be 
understood that the segregation of the various functions into the illustrated blocks 
5 is for the purpose of facilitating a description of the invention. In practice, it is 

possible for some or all of the illustrated blocks to be combined into a single unit 
that performs all of the combined functions. 

In addition to depicting the hardware and structure and flow of 
operations associated with the stage 2 processing (i.e., the fine frame 

10 synchronization), FIG. 26 additionally shows how the iterative synchronization 

procedure of stage 2 may be combined with the overall frequency and frame 
synchronization framework established by the stage 1 processing, which is 
represented within the dashed lines in the figure. That which appears outside the 
dashed lines in FIG. 26 corresponds to stage 2 processing. 

15 The coarse frame start estimate 509 (i.e., k) from a time-domain 

technique 2601 (e.g., according to any of the metrics specified in Eqs. (20) 
through (28)) can be used to extract a frequency offset estimate, g , by determining 
a correlation result, s[k], in accordance with Eq. (18), and then applying Eq. (22) 
to exploit the argument of the correlation result. The frequency offset estimate 

20 would be optimum if the correlation result were evaluated at the correct frame 

start. Due to time offsets, the variance of this frequency estimate will be higher 
than in the optimum case. 

Upon these coarse estimates, the iterative frequency -domain frame 
synchronization unit 2605 starts to work, repeatedly determining a correlation 

25 result and generating therefrom a next time offset estimate | that is used to adjust 

the DFTs demodulation window in the sample memory 501 for use in a next 
iteration until convergence is reached, thereby yielding an optimum frame start 
estimate 2603. In a preferred embodiment, convergence is achieved within two 
iterations, but of course this could be different in other applications. 
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As illustrated in FIG. 26, this optimum frame start estimate 2603 
can again be used to further improve the frequency offset estimate by evaluating 
the previously stored argument at the correct time position within the correlation 
result memory 511. 

5 After convergence of the iterative refinement procedure, both the 

frame synchronization estimate 2603 and the frequency offset estimate 525 are 
reliably obtained with low variance. 

Advantages of the Proposed OFDM Frame Synchronization Scheme 
Advantages associated with the iteratively improved frame 
10 synchronization and the modified OFDM symbol structure in accordance with the 

invention include: 

• The frequency-domain refinement procedure offers a time 
synchronization with very low variance for spontaneous 
transmission over unknown frequency-selective channels, 

15 even though a systematic time synchronization "error" is 

introduced by the channel power delay characteristic. 

• The expected time synchronization position found with this 
technique automatically corrects the phase ramp across the 
subcarriers caused by the center of power delay gravity in 

20 the channel impulse response. This is a very reasonable 

frame synchronization estimate with OFDM, especially in 
conjunction with differential demodulation in the frequency 
direction. 

• The cyclic postfix prevents interference from the following 
25 symbol which would otherwise systematically occur due to 

the channel-delayed frame sync positions. 

• The division of the guard interval into a prefix and a postfix 
offers reduced sensitivity in terms of time synchronization 
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errors in a small interval around the ideal position. More 
concisely, it helps the iterative process to converge more 
easily. 

Simulation Results 

The setup for the presented simulations realizes an OFDM symbol 
structure with a prefix and a postfix guard interval as illustrated in FIG. 25. 
Hence the guard interval of D g = 8 modulation intervals is split into equal parts of 
4 samples representing a cyclic prefix and 4 samples realizing the cyclic postfix. 
Consequently, for transmission over AWGN we will have a misplacement D m = 0 
for sampling offsets |£ t | < 4. The OFDM burst is transmitted over the respective 
channel and the receiver starts the iterative procedure upon a fixed and 
deterministic misplacement D m which represents the parameter in the following 
diagrams. 

In these simulations, the frequency-domain time offset estimation 
was limited to two processing steps, i.e., the iterative frequency domain frame 
synchronization unit 2605 in FIG. 26 was run only twice. Thus, two DFTs were 
required for frame sync and another one for final data demodulation. The 
achieved performance improvement is investigated in terms of the variance of 

the timing estimate % t and in terms of the time synchronization failure rate p < w > 

which is the probability that the estimated frame start position lies more than ±m 
samples away from the actual optimum frame start. For OFDM transmission, 
depending on the number of subcarriers £>, it can be sufficient to choose anm^O. 
In a system having 64 carriers, we choose m = 1, so that there is an interval of 3 
timing positions which are interpreted as valid (acceptable) symbol starts. Hence 
we allow, at most, one interfering sample in the demodulation window of length D 
= 64 resulting in minimally 18.06 dB SNR due to frame misalignment. 
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For the simulation of the estimation variance, the estimator 
parameters D cg and D c - lg present in the lower bound of Eq. (70) were varied. The 
variation was done in such a way so as to constantly maintain D cg - D cig = 10, 
i.e., a constant training symbol overhead was maintained in order to enable a fair 
5 comparison. The simulated standard deviations for two parameter settings and 

transmission over an AWGN channel are depicted in FIGS. 27 and 28. More 
specifically, FIGS. 27 and 28 are each graphs showing the simulated standard 
deviation of £ in the iterative frame synchronization for transmission over AWGN 
in the presence of an initial misplacement D m . In each case, the lower bound is 
10 according to Eq. (70). In the graph of FIG. 27, the frequency-domain estimator 

worked with £> cg = 10 and D cig = 1 in a system with D= 64, £> g =8 and £> sync =32. 
In the graph of FIG. 28, the frequency-domain estimator worked with D cg =5 and " 
£> cig =2 in a system with £>=64, D g =8 and Z> sync =32. 

The first thing to notice is the constantly decreasing lower bound for 
15 the time offset estimate when the number of correlations in one group D cig is 

increased. This behavior can analytically be seen from Eq. (70) and it is due to 
the summation of adjacent frequency-domain correlation products in the 
overlapping estimators of the type £ 0 = 1 . It can be seen from the graphs of 
FIGS. 27 and 28 that up to a misplacement of D m = 4 (i.e., | § t | < 8), the 
20 simulated estimation variance gets increasingly lower. This means that for 

moderate time offsets, the simulation results manage to follow the decreasing 
lower bound, even though it becomes visibly more problematic, i.e., for all D m , 
the asymptotic behavior is not reached until higher SNRs. 

The situation is quite different for larger initial time offsets. 
25 Generally there is an improvement for the step from FIG. 27 (D cg = 10 and D cig 

= 1) to FIG. 28 (D cg = 5 and D cig = 2). However, a further step to (D cg = 2 
and Z) cig = 5) would result in a degradation for the (iterated) estimation 
performance for D m ^ 6. In conclusion, the variant D cg = 5 and £> cig = 2 
embodies the most robust of the investigated parameter settings. 
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Turning now to the graphs of FIGS. 29 and 30, we consider the 
time synchronization failure rates for the previously investigated AWGN 
transmission scenario, that is, transmission over AWGN in the presence of an 
initial misplacement D m . In the graph of FIG. 29, the frequency-domain estimator 
5 works with D cg = 10 and D cig =l in a system with £>=64, D g =8 and D sync = 32. 

In the graph of FIG. 30, the frequency-domain estimator works with D cg =5 and 
D cig =2 in a system with D=64, D g = 8 and jD sync =32. Here we observe the same 
estimator behavior as was noted above with respect to FIGS. 27 and 28: Moving 
from the estimator parameters in FIG. 29 to the ones in FIG. 30 results in an 

10 overall improvement of the performance. Gains of 2 dB and more are realized. 

The AWGN performance in the latter diagram reveals that two iterations of the 
iterative frequency-domain approach are able to transform an initial misplacement 
of D m = 8 (i.e., |£ t | < 12) into a successful time synchronization (i.e., £ -£) e 
{-1,0, -f 1}) in over 99% of cases for 10 log 10 (E s /N 0 ) ;> 6 dB. Having in mind 

15 that this estimator operates with only £> sync = 32, the results are quite impressive. 

Larger OFDM symbols will result in larger allowable initial offsets, i.e., a larger 
range of convergence and more iterations would result in lower error floors for 
p( m ). 

Similar performance gains were also obtained for transmission over 
20 a multipath channel, where D cg = 5 and D cig = 2 was found to be a good choice. 

In frequency-selective channels, D cg must necessarily be quite large, as it 
determines the diversity factor achieved by the frequency-domain pilot symbols, 
and thus it directly affects the power efficiency of the iterative frame 
synchronization performance. 
25 One can conclude from the simulation results for the iterative 

frequency-domain frame synchronization applied to the exemplary system that two 
iterated frequency-domain estimates are adequate for providing a robust frame 
synchronization refinement for initial coarse time estimates in a range of at least 
|£ t | < 6 for multipath and |£ t | < 10 for AWGN channels. For moderate offsets, 
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the remaining error floor in the synchronization failure rate diagrams can be 
further reduced by an increased number of iterations. For very large offsets as 
starting points of the iteration, the error floor is caused by the first processing loop 
which eventually drives the frame entirely out of the region of convergence. 

5 OVERALL SYNCHRONIZATION PERFORMANCE 

It is desired to evaluate the Packet Error Rate (PER) over SNR in a 

severe multipath channel with operational synchronization in a D — 64 carrier 

system. The synchronization has been chosen to work with a sandamble with 

£>«mn = 32. If not otherwise stated, the historic MMSE criterion according to Eq. 
sync 

10 (21) is used for coarse frame synchronization, even though the MNC metric of Eq. 

(28) was found to be more robust and preferable. 

For frame synchronization refinement, only the first frequency- 
domain frame synchronization estimate is used, and the results are not iterated. 
Here, we use Z? cg = 5 correlation groups with D cig = 2 correlations in one group, 

15 i.e., 10 dedicated differential pilots are introduced into the frequency multiplex of 

the sandamble synchronization symbol. 

For the simulations a constant NCFO of £ f = 0.3 is used. This is 
beyond the estimation range of the pure sandamble correlation. Hence, a FD 
frequency detector with D c = 10 and 3 discrete intervals is used that operates with 

20 uniform a priori probabilities for all three detected frequency offset intervals. 

FIG. 31 is a graph depicting, for postfix symbols comprising 
various numbers of samples (i.e., 0, 1, 2 and 3 samples per postfix symbol), the 
PER as a function of SNR for convolutionally coded 8-Differential Phase Shift 
Keying (8DPSK) modulation in the time direction for ideal synchronization and 

25 operational burst synchronization before and after the frame synchronization 

refinement. The PER represents the performance after Soft Decision (SD) Viterbi- 
Decoding. The packet error rate indicates the rate of received packets with 
residual transmission errors. Each packet consists of 570 user bits. From this 
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graph one concludes that performance improves when a postfix of reasonable 
length is used. That is, the use of no postfix as well as the use of too much postfix 
yields poor results. Using a postfix of 2 samples is found to be a good 
compromise in this severe multipath channel. The performance achieves the ideal 
5 PER within a fraction of 1 dB. 

FIG. 32 provides an idea of the consequences of an unreliable 
coarse timing estimate. More specifically, FIG. 32 is a graph of the PERs plotted 
as a function of SNR, where PERs are determined for the cases of ideal 
synchronization and operational burst synchronization before and after the frame 

10 synchronization refinement. The presented results are for 8DPSK modulation in 

the time direction. However, similar results are obtained for the frequency 
direction as well. Here, the frequency-domain frame sync refinement shows large 
gains, but it cannot repair the major problems caused by the bad initial estimate. 
The bad initial estimate from the MC criterion ruins the overall performance; the 

15 PER after the frame sync refinement hits an error floor around 10 . 

FIG. 33 is a graph showing a comparison between differential 
modulation in the frequency direction and in the time direction. For each, PERs 
are plotted as a function of SNR for ideal synchronization and operational burst 
synchronization before and after the frame synchronization refinement. In this 

20 diagram there are three ideal PERs (graphs 3301, 3303 and 3305), with one of 

them (graph 3305) being for 8DPSK modulation in the time direction. The other 
two ideal curves (graphs 3301 and 3303) are for 8DPSK modulation in the 
frequency direction. The lower of these (graph 3303) is with an MMSE optimal 
demodulation of differential modulation schemes in the frequency direction. The 

25 higher of these (graph 3301) is for "ideal" demodulation window alignment with 

no channel phase correction. 

Two important conclusions can be drawn from FIG. 33: 
First, for DPSK in the frequency direction with small numbers of 
carriers (here D = 64), a frequency domain frame sync refinement is called for 
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because the performance with the coarse time estimate (indicated by graph 3307) is 
not satisfying. Furthermore, the performance with refinement (indicated by graph 
3309) approaches that of the MMSE demodulation (graph 3303) within a fraction 
of one dB. This confirms the derivation presented earlier. 

Also, for DPSK in the time direction, the performance with the 



coarse time estimate (indicated by graph 3311) is already satisfying. As indicated 

_3 

by graph 3313, a performance gain of 1 dB at PER = 10 is possible with the 
frequency-domain frame sync refinement. 



SUMMARY OF FEATURES 
10 A number of aspects of the invention will be apparent to those 

having ordinary skill in the art. These include: 
In the context of frequency synchronization: 

The sandamble structure itself. This includes the temporal division 
of the periodically repeated signal parts in two (or more) regions, 
15 eventually separated by data. 

The use of an extended guard region within a sandamble to make it 
possible to resolve frequency ambiguities by the FD approach. 

The FD approach for frequency offset estimation. This implies the 
first (fine) estimation due to the correlation between A and B and 
20 the correction of this estimated offset by techniques such as those 

illustrated in FIGS. 5a and 5b. 

A subsequent correlation between G 2 and A 2 and/or G 2 and B 2 (and 
in at least one embodiment, correction of these correlation results) and the 
advantageous combination (e.g., addition) of these two (possibly corrected) 
25 results. Detection (decision) of the residual frequency offset interval based on the 

final correlation result (e.g., combined correlation result). 
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The use of a priori knowledge (probabilities) of the possible 
frequency offsets in the detection step of FD to optimize the overall 
frequency synchronization performance. 

In the context of frame (time) synchronization: 
5 - The guard interval, which usually consists of only a cyclic prefix, is 

supplemented with a cyclic postfix. 

The cyclic prefix can be shortened, and this spare time added to the 
end of the useful OFDM symbol part for use as the cyclic postfix. 

In the system design, the distribution of the total number of guard 

10 samples in the pre- and postfix can be optimally adapted to the expected channel 

power delay profile. 

The frame synchronization refinement can start upon a coarse 
estimate. The only requirement is that it must be reliable enough to 
ensure convergence of the frequency-domain based estimator. 

15 - The estimate for frame synchronization may be obtained from a 

frequency-domain processing stage, which operates iteratively. 
Each iteration step includes 1) frame alignment according to the 
most recent estimate, 2) either a complete DFT or a recursive 
sliding DFT of this frame, 3) time offset estimation from the 

20 frequency-domain amplitudes. If the algorithm converges, the 

variance of this estimate will get increasingly smaller. 
Of course, only one single step of the proposed frequency-domain 
refinement need be implemented (and therefore executed), so that in 
this case the procedure is non-iterative. In other words, one is free 

25 to perform any number of refinements (including only one 

refinement), as appropriate. 

The iterative processing stage need not operate on a special 
preamble symbol. In a decision-directed approach it can be applied 
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to every OFDM symbol in the burst. The pilot symbols in a data- 
aided approach can be multiplexed into any OFDM symbol in the 
burst. When the frequency offset is already corrected and a reliable 
coarse timing estimate can be derived from periodic protocol 
structures in the radio link, the second stage can be used in a stand- 
alone mode without the sandamble. In this case we arrive at a very 
efficient "preambieless" frame synchronization scheme. 
In a combination of stages 1 and 2, the symbols in a data-aided 
approach can either be multiplexed into any OFDM symbol in the 
burst or into the sandamble symbol. Depending on the specific 
system parameters one of these approaches may be advantageous 
than the other. 

The last frame synchronization estimate (or even that after each 
iteration) can be used to steadily improve the frequency offset 
correction. After the described synchronization procedure 
performed upon the burst preamble or the first symbol, the 
remaining OFDM burst can be demodulated with nearly optimum 
(minimum variance) frequency and time offset estimates. 

The invention has been described with reference to a particular 
embodiment. However, it will be readily apparent to those skilled in the art that it 
is possible to embody the invention in specific forms other than those of the 
preferred embodiment described above. This may be done without departing from 
the spirit of the invention. The preferred embodiment is merely illustrative and 
should not be considered restrictive in any way. The scope of the invention is 
given by the appended claims, rather than the preceding description, and all 
variations and equivalents which fall within the range of the claims are intended to 
be embraced therein. 
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WHAT IS CLAIMED IS: 

1 . A method of transmitting a sequence of data samples, comprising 
the steps of: 

initially transmitting a first preamble comprising a sequence of 
5 preamble samples; 

transmitting the sequence of data samples; and 
subsequently transmitting a second preamble comprising the 
sequence of preamble samples; 

whereby the sequence of data samples is transmitted after the step of 
10 initially transmitting the first preamble, and before the step of subsequently 

transmitting the second preamble. 

2. The method of claim 1, wherein: 

the first preamble comprises a first guard region and a symbol 

sequence; 

15 the first guard region comprises a first subset of the symbol 

sequence; 

the second preamble comprises a second guard region and the 
symbol sequence; and 

the second guard region comprises a second subset of the symbol 

20 sequence. 



3. The method of claim 2, wherein the first subset of the symbol 

sequence is a larger subset of the symbol sequence than the second subset of the 
symbol sequence. 
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4. The method of claim 2, wherein the second subset of the symbol 
sequence is a larger subset of the symbol sequence than the first subset of the 
symbol sequence. 

5. A method of receiving a sequence of desired data samples, 
5 comprising the steps of: 

receiving a signal sequence that includes a first preamble followed 
by the sequence of desired data samples followed by a second preamble, wherein 
the first preamble comprises a symbol sequence, and the second preamble 
comprises the symbol sequence; 
10 hypothesizing a plurality of locations of the first preamble in the 

received plurality of data samples; 

for each of the hypothesized locations, determining a hypothesized 
first preamble and a corresponding hypothesized second preamble; 

for each of the hypothesized locations, determining a correlation 
1 5 between the hypothesized first preamble and the corresponding hypothesized 

second preamble; and 

using the correlations to determine locations of a most likely 
hypothesized first preamble and of a most likely second preamble. 

6. The method of claim 5, further comprising the step of: 

20 using the locations of the most likely hypothesized first and second 

preambles to estimate a first frequency offset of the received signal sequence. 

7. The method of claim 6, wherein: 

the first preamble comprises a first guard region and a symbol 

sequence; 

25 the first guard region comprises a first subset of the symbol 

sequence and a second subset of the symbol sequence; 
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the second preamble comprises a second guard region and the 
symbol sequence; and 

the method further comprises the steps of: 

using the first frequency offset to make a first frequency 
5 correction of the received signal sequence; 

determining a frequency corrected first guard region in the 
first frequency corrected received signal sequence; 

determining a frequency corrected first preamble in the first 
frequency corrected received signal sequence; 
10 correlating the second subset of the symbol sequence from 

the frequency corrected first guard region with the second subset of the symbol 
sequence from the frequency corrected first preamble, thereby generating a 
correlation result; and 

using the correlation result to detect a second frequency 

15 offset. 



8. The method of claim 7, wherein the step of using the correlation 

result to detect the second frequency offset comprises the step of using a priori 
knowledge of probabilities of possible frequency offset intervals to detect the 
second frequency offset. 

20 9. The method of claim 6, wherein: 

the first preamble comprises a first guard region and a symbol 

sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 
25 the second preamble comprises a second guard region and the 

symbol sequence; and 

the method further comprises the steps of: 
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correlating the second subset of the symbol sequence from 
the first guard region with the second subset of the symbol sequence from the first 
preamble, thereby generating a correlation result; 

using the first frequency offset to correct the correlation 

5 result; and 

using the corrected correlation result to estimate a second 

frequency offset. 

10. The method of claim 6, wherein: 

the first preamble comprises a first guard region and a symbol 

10 sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 

the second preamble comprises a second guard region and the 
symbol sequence; and 
15 the method further comprises the steps of: 

using the first frequency offset to make a first frequency 
correction of the received signal sequence; 

determining a frequency corrected first guard region in the 
first frequency corrected received signal sequence; 
20 determining a frequency corrected second preamble in the 

first frequency corrected received signal sequence; 

correlating the second subset of the symbol sequence from 
the frequency corrected first guard region with the second subset of the symbol 
sequence from the frequency corrected second preamble, thereby generating a 
25 correlation result; and 

using the correlation result to detect a second frequency 

offset. 
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11. The method of claim 10, wherein the step of using the correlation 

result to detect the second frequency offset comprises the step of using a priori 
knowledge of probabilities of possible frequency offset intervals to detect the 
second frequency offset. 

5 12. The method of claim 6, wherein: 

the first preamble comprises a first guard region and a symbol 

sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 
10 the second preamble comprises a second guard region and the 

symbol sequence; and 

the method further comprises the steps of: 

correlating the second subset of the symbol sequence from 
the first guard region with the second subset of the symbol sequence from the 
15 second preamble, thereby generating a correlation result; 

using the first frequency offset to correct the correlation 

result; and 

using the corrected correlation result to estimate a second 

frequency offset. 



20 13. The method of claim 6, wherein: 

the first preamble comprises a first guard region and a symbol 

sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 
25 the second preamble comprises a second guard region and the 

symbol sequence; and 

the method further comprises the steps of: 
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using the first frequency offset to make a first frequency 
correction of the received signal sequence; 

determining a frequency corrected first guard region in the 
first frequency corrected received signal sequence; 
5 determining a frequency corrected first preamble in the first 

frequency corrected received signal sequence; 

determining a corrected second preamble in the first 
frequency corrected received signal sequence; 

correlating the second subset of the symbol sequence from 
10 the frequency corrected first guard region with the second subset of the symbol 

sequence from the frequency corrected first preamble, thereby generating a first 
correlation result; 

correlating the second subset of the symbol sequence from 
the frequency corrected first guard region with the second subset of the symbol 
15 sequence from the frequency corrected second preamble, thereby generating a 

second correlation result; 

combining the first and second correlation results, thereby 
generating a combined correlation result; and 

using the combined correlation result to detect a second 

20 frequency offset. 

14. The method of claim 13, wherein the step of using the combined 

correlation result to detect the second frequency offset comprises the step of using 
a priori knowledge of probabilities of possible frequency offset intervals to detect 
the second frequency offset. 

25 15. The method of claim 6, wherein: 

the first preamble comprises a first guard region and a symbol 



sequence; 
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the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 

the second preamble comprises a second guard region and the 
symbol sequence; and 
5 the method further comprises the steps of: 

correlating the second subset of the symbol sequence from 
the first guard region with the second subset of the symbol sequence from the first 
preamble, thereby generating a first correlation result; 

correlating the second subset of the symbol sequence from 
10 the first guard region with the second subset of the symbol sequence from the 

second preamble, thereby generating a second correlation result; 

using the first frequency offset to correct the first correlation 

result; 

using the first frequency offset to correct the second 

1 5 correlation result ; 

combining the first and second corrected correlation results, 
thereby generating a combined corrected correlation result; and 

using the combined corrected correlation result to estimate a 
second frequency offset. 

20 16. The method of claim 5, wherein: 

the received signal sequence includes the first preamble, followed 
by a cyclic prefix guard interval, followed by the sequence of desired data 
samples, followed by a cyclic postfix guard interval, followed by a second 
preamble, and 

25 the method further comprises the steps of: 

using the location of the most likely hypothesized first preamble to 
determine a coarse estimate of a starting position of the sequence of desired 
samples in the received signal sequence; and 
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generating an improved estimate of the starting position of the 
sequence of desired samples in the received signal sequence by: 

using a time domain to frequency domain transformation 
technique and the coarse estimate of the starting position of the sequence of desired 
5 samples in the received signal sequence to generate frequency-domain received 

samples; 

determining a frequency-domain correlation between the 
frequency-domain received samples and noiseless samples; and 

generating the improved estimate of the starting position of 
10 the sequence of desired samples in the received signal sequence based on an 

argument of the frequency-domain correlation between the frequency-domain 
received samples and the noiseless samples. 

17. The method of claim 16, wherein: 

the cyclic prefix guard interval comprises an ending subset of the 
15 sequence of desired data samples; and 

the cyclic postfix guard interval comprises a beginning subset of the 
sequence of desired data samples. 

18. The method of claim 16, wherein the number of samples in the 
cyclic prefix guard interval and the number of samples in the cyclic postfix guard 

20 interval are optimally adapted to an expected channel power delay profile. 

19. The method of claim 16, further comprising the steps of: 
generating a further improved estimate of the starting position of the 

sequence of desired samples in the received signal sequence by: 

using the improved estimate of the starting position of the 
25 sequence of desired samples in the received signal sequence to generate improved 

frequency -domain received samples; 
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detennining a second frequency-domain correlation between 
the improved frequency-domain received samples and the noiseless samples; and 

generating the further improved estimate of the starting 
position of the sequence of desired samples in the received signal sequence based 
5 on an argument of the second frequency-domain correlation between the 

frequency-domain received samples and the noiseless samples. 

20. The method of claim 16, further comprising the step of generating 

the noiseless samples by demodulating symbols contained in the received signal 
sequence. 



10 21. The method of claim 20, further comprising the step of generating 

the noiseless samples by demodulating symbols contained in the sequence of 
desired samples. 

22. The method of claim 16, further comprising the step of generating 

the noiseless samples from pilot symbols that have been multiplexed into the 
15 sequence of received samples. 



23. A method of receiving a sequence of desired data samples, 

comprising the steps of: 

receiving a signal sequence that includes a cyclic prefix guard 
interval followed by the sequence of desired data samples, followed by a cyclic 
20 postfix guard interval; 

determining a coarse estimate of a starting position of the sequence 
of desired samples in the received signal sequence; and 

generating an improved estimate of the starting position of the 
sequence of desired samples in the received signal sequence by: 
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using the coarse estimate of the starting position of the 
sequence of desired samples in the received signal sequence to generate frequency- 
domain received samples; 

determining a frequency-domain correlation between the 
5 frequency-domain received samples and noiseless samples; and 

generating the improved estimate of the starting position of 
the sequence of desired samples in the received signal sequence based on an 
argument of the frequency -domain correlation between the frequency-domain 
received samples and the noiseless samples. 

10 24. The method of claim 23, wherein: 

the cyclic prefix guard interval comprises an ending subset of the 
sequence of desired data samples; and 

the cyclic postfix guard interval comprises a beginning subset of the 
sequence of desired data samples. 

15 25. The method of claim 23, wherein the number of samples in the 

cyclic prefix guard interval and the number of samples in the cyclic postfix guard 
interval are optimally adapted to an expected channel power delay profile. 

26. The method of claim 23, further comprising the steps of: 

generating a further improved estimate of the starting position of the 
20 sequence of desired samples in the received signal sequence by: 

using the improved estimate of the starting position of the 
sequence of desired samples in the received signal sequence to generate improved 
frequency -domain received samples; 

determining a second frequency-domain correlation between 
25 the improved frequency -domain received samples and the noiseless samples; and 
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generating the further improved estimate of the starting 
position of the sequence of desired samples in the received signal sequence based 
on an argument of the second frequency-domain correlation between the 
frequency-domain received samples and the noiseless samples. 



5 27. The method of claim 26, further comprising the step of using the 

further improved estimate of the starting position of the sequence of desired 
samples in the received signal sequence to adjust a measurement of frequency 
offset associated with the sequence of desired samples in the received signal 
sequence. 



10 28. The method of claim 23, wherein the noiseless samples are 

generated by demodulating symbols contained in the received signal sequence. 

29. The method of claim 28, wherein the noiseless samples are 
generated by demodulating symbols contained in the sequence of desired samples. 

30. The method of claim 23, wherein the noiseless samples are 

15 generated from pilot symbols that have been multiplexed into the sequence of 

received samples. 



31. The method of claim 23, further comprising the step of using the 

improved estimate of the starting position of the sequence of desired samples in the 
received signal sequence to adjust a measurement of frequency offset associated 
20 with the sequence of desired samples in the received signal sequence. 



32. The method of claim 23, further comprising the steps of: 

iteratively generating a final improved estimate of the starting 
position of the sequence of desired samples in the received signal sequence by 



WO 99/65180 



PCT/EP99/03793 



-99- 

performing a sequence of iterative steps for a number of times, each sequence of 
iterative steps comprising: 

using a previously generated estimate of the starting position 
of the sequence of desired samples in the received signal sequence to generate 
5 current frequency-domain received samples; 

determining a current frequency-domain correlation between 
the current frequency-domain received samples and the noiseless samples; and 

generating a next estimate of the starting position of the 
sequence of desired samples in the received signal sequence based on an argument 
10 of the current frequency-domain correlation between the frequency -domain 

received samples and the noiseless samples. 

33. The method of claim 32, further comprising the step of using the 
final improved estimate of the starting position of the sequence of desired samples 
in the received signal sequence to adjust a measurement of frequency offset 

15 associated with the sequence of desired samples in the received signal sequence. 

34. An apparatus for transmitting a sequence of data samples, 
comprising: 

means for initially transmitting a first preamble comprising a 
sequence of preamble samples; 
20 means for transmitting the sequence of data samples; and 

means for subsequently transmitting a second preamble comprising 
the sequence of preamble samples; 

whereby the sequence of data samples is transmitted after the step of 
initially transmitting the first preamble, and before the step of subsequently 
25 transmitting the second preamble. 



35. 



The apparatus of claim 34, wherein: 
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the first preamble comprises a first guard region and a symbol 

sequence; 

the first guard region comprises a first subset of the symbol 

sequence; 

5 the second preamble comprises a second guard region and the 

symbol sequence; and 

the second guard region comprises a second subset of the symbol 

sequence. 

36. The apparatus of claim 35, wherein the first subset of the symbol 
10 sequence is a larger subset of the symbol sequence than the second subset of the 

symbol sequence. 

37. The apparatus of claim 35, wherein the second subset of the symbol 
sequence is a larger subset of the symbol sequence than the first subset of the 
symbol sequence. 

15 38. An apparatus for receiving a sequence of desired data samples, 

comprising: 

means for receiving a signal sequence that includes a first preamble 
followed by the sequence of desired data samples followed by a second preamble, 
wherein the first preamble comprises a symbol sequence, and the second preamble 
20 comprises the symbol sequence; 

means for hypothesizing a plurality of locations of the first preamble 
in the received plurality of data samples; 

means for determining a hypothesized first preamble and a 
corresponding hypothesized second preamble for each of the hypothesized 
25 locations; 
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means for determining a correlation between the hypothesized first 
preamble and the corresponding hypothesized second preamble for each of the 
hypothesized locations; and 

means for using the correlations to determine locations of a most 
5 likely hypothesized first preamble and of a most likely second preamble. 

39. The apparatus of claim 38, further comprising: 

means for using the locations of the most likely hypothesized first 
and second preambles to estimate a first frequency offset of the received signal 
sequence. 

10 40. The apparatus of claim 39, wherein: 

the first preamble comprises a first guard region and a symbol 

sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 
15 the second preamble comprises a second guard region and the 

symbol sequence; and 

the apparatus further comprises: 

means for using the first frequency offset to make a first 
frequency correction of the received signal sequence; 
20 means for determining a frequency corrected first guard 

region in the first frequency corrected received signal sequence; 

means for determining a frequency corrected first preamble 
in the first frequency corrected received signal sequence; 

means for correlating the second subset of the symbol 
25 sequence from the frequency corrected first guard region with the second subset of 

the symbol sequence from the frequency corrected first preamble, thereby 
generating a correlation result; and 
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means for using the correlation result to detect a second 

frequency offset. 

41. The apparatus of claim 40, wherein the means for using the 
correlation result to detect the second frequency offset comprises means for using a 

5 priori knowledge of probabilities of possible frequency offset intervals to detect the 

second frequency offset. 

42. The apparatus of claim 39, wherein: 
the first preamble comprises a first guard region and a symbol 

sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 

the second preamble comprises a second guard region and the 
symbol sequence; and 

the apparatus further comprises: 

means for correlating the second subset of the symbol 
sequence from the first guard region with the second subset of the symbol 
sequence from the first preamble, thereby generating a correlation result; 

means for using the first frequency offset to correct the 
correlation result; and 

means for using the corrected correlation result to estimate a 
second frequency offset. 

43. The apparatus of claim 39, wherein: 
the first preamble comprises a first guard region and a symbol 

sequence; 

25 the first guard region comprises a first subset of the symbol 

sequence and a second subset of the symbol sequence; 



10 



15 



20 
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the second preamble comprises a second guard region and the 
symbol sequence; and 

the apparatus further comprises: 

means for using the first frequency offset to make a first 
5 frequency correction of the received signal sequence; 

means for determining a frequency corrected first guard 
region in the first frequency corrected received signal sequence; 

means for determining a frequency corrected second 
preamble in the first frequency corrected received signal sequence; 
10 means for correlating the second subset of the symbol 

sequence from the frequency corrected first guard region with the second subset of 
the symbol sequence from the frequency corrected second preamble, thereby 
generating a correlation result; and 

means for using the correlation result to detect a second 

15 frequency offset. 

44. The apparatus of claim 43, wherein the means for using the 

correlation result to detect the second frequency offset comprises means for using a 
priori knowledge of probabilities of possible frequency offset intervals to detect the 
second frequency offset. 

20 45. The apparatus of claim 39, wherein: 

the first preamble comprises a first guard region and a symbol 

sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 
25 the second preamble comprises a second guard region and the 

symbol sequence; and 

the apparatus farther comprises: 
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means for correlating the second subset of the symbol 
sequence from the first guard region with the second subset of the symbol 
sequence from the second preamble, thereby generating a correlation result; 

means for using the first frequency offset to correct the 

5 correlation result; and 

means for using the corrected correlation result to estimate a 
second frequency offset. 

46. The apparatus of claim 39, wherein: 

the first preamble comprises a first guard region and a symbol 

10 sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 

the second preamble comprises a second guard region and the 
symbol sequence; and 
15 the apparatus further comprises: 

means for using the first frequency offset to make a first 
frequency correction of the received signal sequence; 

means for determining a frequency corrected first guard 
region in the first frequency corrected received signal sequence; 
20 means for determining a frequency corrected first preamble 

in the first frequency corrected received signal sequence; 

means for determining a corrected second preamble in the 
first frequency corrected received signal sequence; 

means for correlating the second subset of the symbol 
25 sequence from the frequency corrected first guard region with the second subset of 

the symbol sequence from the frequency corrected first preamble, thereby 
generating a first correlation result; 
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means for correlating the second subset of the symbol 
sequence from the frequency corrected first guard region with the second subset of 
the symbol sequence from the frequency corrected second preamble, thereby 
generating a second correlation result; 
5 means for combining the first and second correlation results, 

thereby generating a combined correlation result; and 

means for using the combined correlation result to detect a 
second frequency offset. 

47. The apparatus of claim 46, wherein the means for using the 

10 combined correlation result to detect the second frequency offset comprises means 

for using a priori knowledge of probabilities of possible frequency offset intervals 
to detect the second frequency offset. 



48. The apparatus of claim 39, wherein: 

the first preamble comprises a first guard region and a symbol 

15 sequence; 

the first guard region comprises a first subset of the symbol 
sequence and a second subset of the symbol sequence; 

; the second preamble comprises a second guard region and the 
symbol sequence; and 
20 the apparatus further comprises: 

means for correlating the second subset of the symbol 
sequence from the first guard region with the second subset of the symbol 
sequence from the first preamble, thereby generating a first correlation result; 

means for correlating the second subset of the symbol 
25 sequence from the first guard region with the second subset of the symbol 

sequence from the second preamble, thereby generating a second correlation 
result; 
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means for using the first frequency offset to correct the first 

correlation result; 

means for using the first frequency offset to correct the 
second correlation result; 
5 means for combining the first and second corrected 

correlation results, thereby generating a combined corrected correlation result; and 

means for using the combined corrected correlation result to 
estimate a second frequency offset. 

49. The apparatus of claim 38, wherein: 

10 the received signal sequence includes the first preamble, followed 

by a cyclic prefix guard interval, followed by the sequence of desired data 
samples, followed by a cyclic postfix guard interval, followed by a second 
preamble, and 

the apparatus further comprises: 
15 means for using the location of the most likely hypothesized first 

preamble to determine a coarse estimate of a starting position of the sequence of 

desired samples in the received signal sequence; and 

means for generating an improved estimate of the starting position 

of the sequence of desired samples in the received signal sequence by: 
20 using a time domain to frequency domain transformation 

technique and the coarse estimate of the starting position of the sequence of desired 

samples in the received signal sequence to generate frequency-domain received 

samples; 

determining a frequency-domain correlation between the 
25 frequency-domain received samples and noiseless samples; and 

generating the improved estimate of the starting position of 
the sequence of desired samples in the received signal sequence based on an 
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argument of the frequency-domain correlation between the frequency -domain 
received samples and the noiseless samples. 

50. The apparatus of claim 49, wherein: 

the cyclic prefix guard interval comprises an ending subset of the 
sequence of desired data samples; and 

the cyclic postfix guard interval comprises a beginning subset of the 
sequence of desired data samples. 

51. The apparatus of claim 49, wherein the number of samples in the 
cyclic prefix guard interval and the number of samples in the cyclic postfix guard 
interval are optimally adapted to an expected channel power delay profile. 

52. The apparatus of claim 49, further comprising: 

means for generating a further improved estimate of the starting 
position of the sequence of desired samples in the received signal sequence by: 

using the improved estimate of the starting position of the 
sequence of desired samples in the received signal sequence to generate improved 
frequency-domain received samples; 

determining a second frequency-domain correlation between 
the improved frequency-domain received samples and the noiseless samples; and 

generating the further improved estimate of the starting 
position of the sequence of desired samples in the received signal sequence based 
on an argument of the second frequency-domain correlation between the 
frequency-domain received samples and the noiseless samples. 

53. The apparatus of claim 49, further comprising means for generating 
the noiseless samples by demodulating symbols contained in the received signal 
sequence. 



WO 99/65180 



PCT/EP99/03793 



-108- 

54. The apparatus of claim 53, further comprising means for generating 
the noiseless samples by demodulating symbols contained in the sequence of 
desired samples. 

55. The apparatus of claim 49, further comprising means for generating 
the noiseless samples from pilot symbols that have been multiplexed into the 
sequence of received samples. 

56. An apparatus for receiving a sequence of desired data samples, 
comprising: 

means for receiving a signal sequence that includes a cyclic prefix 
guard interval followed by the sequence of desired data samples, followed by a 
cyclic postfix guard interval; 

means for determining a coarse estimate of a starting position of the 
sequence of desired samples in the received signal sequence; and 

means for generating an improved estimate of the starting position 
of the sequence of desired samples in the received signal sequence by: 

using the coarse estimate of the starting position of the 
sequence of desired samples in the received signal sequence to generate frequency- 
domain received samples; 

determining a frequency-domain correlation between the 
frequency-domain received samples and noiseless samples; and 

generating the improved estimate of the starting position of 
the sequence of desired samples in the received signal sequence based on an 
argument of the frequency-domain correlation between the frequency-domain 
received samples and the noiseless samples. 



57. 



The apparatus of claim 56, wherein: 
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the cyclic prefix guard interval comprises an ending subset of the 
sequence of desired data samples; and 

the cyclic postfix guard interval comprises a beginning subset of the 
sequence of desired data samples. 

5 58. The apparatus of claim 56, wherein the number of samples in the 

cyclic prefix guard interval and the number of samples in the cyclic postfix guard 
interval are optimally adapted to an expected channel power delay profile. 

59. The apparatus of claim 56, further comprising: 

means for generating a further improved estimate of the starting 
10 position of the sequence of desired samples in the received signal sequence by: 

using the improved estimate of the starting position of the 
sequence of desired samples in the received signal sequence to generate improved 
frequency-domain received samples; 

determining a second frequency-domain correlation between 
15 the improved frequency-domain received samples and the noiseless samples; and 

generating the further improved estimate of the starting 
position of the sequence of desired samples in the received signal sequence based 
on an argument of the second frequency-domain correlation between the 
frequency-domain received samples and the noiseless samples. 

20 60. The apparatus of claim 59, further comprising means for using the 

further improved estimate of the starting position of the sequence of desired 
samples in the received signal sequence to adjust a measurement of frequency 
offset associated with the sequence of desired samples in the received signal 
sequence. 
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61. The apparatus of claim 56, further comprising means for generating 
the noiseless samples by demodulating symbols contained in the received signal 
sequence. 

62. The apparatus of claim 61, further comprising means for generating 
the noiseless samples by demodulating symbols contained in the sequence of 
desired samples. 

63. The apparatus of claim 56, further comprising means for generating 
the noiseless samples from pilot symbols that have been multiplexed into the 
sequence of received samples. 

64. The apparatus of claim 56, further comprising means for using the 
improved estimate of the starting position of the sequence of desired samples in the 
received signal sequence to adjust a measurement of frequency offset associated 
with the sequence of desired samples in the received signal sequence. 

65. The apparatus of claim 56, further comprising: 

means for iteratively generating a final improved estimate of the 
starting position of the sequence of desired samples in the received signal sequence 
by performing a sequence of iterative steps for a number of times, each sequence 
of iterative steps comprising: 

using a previously generated estimate of the starting position 
of the sequence of desired samples in the received signal sequence to generate 
current frequency -domain received samples; 

determining a current frequency-domain correlation between 
the current frequency -domain received samples and the noiseless samples; and 

generating a next estimate of the starting position of the 
sequence of desired samples in the received signal sequence based on an argument 
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of the current frequency-domain correlation between the frequency-domain 
received samples and the noiseless samples. 

66. The apparatus of claim 65, further comprising means for using the 

final improved estimate of the starting position of the sequence of desired samples 
5 in the received signal sequence to adjust a measurement of frequency offset 

associated with the sequence of desired samples in the received signal sequence. 



WO 99/65180 



1/31 



PCT/EP99/03793 



FIG. 1a 

(PRIOR ART) 



G 


A 

A 


B 


data 




0 101 


FIG. 1b 

(PRIOR ART) 


G 


A 


B 


C 




s 

103 







data 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



2/31 



PCT/EP99/03793 




SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



PCT/EP99/03793 




SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



PCT/EP99/03793 



4/31 



as 
CD 



CSJ 
GO 



CO 
CD 



CO 
CD 



CXI 



CNJ 
CD 

cd" 



CO 
LO 

L 



CD ' 
LO 



L -T 

o 
E 

CD 

E 

CD 



ca 

CO 



o 

CD 



o 
o 



o 
cz 

CD 



CD 



CO 



CD 
LO 



.CZ5 
LO 



LO 
CD ' 
LO 





ca 




CO 


meti 


CO 
CD 
"O 


"ca 


cn 

CO 




=3 


CD 


CD 


O 






E 


o 








o 





LO 



a 







.2 




"ca 




CD 


§1 

CD CO 


o 


E <c 


CJ> 






SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



PCT/EP99/03793 



5/31 



LO 

CD 



CNJ | 
CO 



QQ 



CO 
CD 



CNJ 
CD 



o 
E 

CD 



o> 

CZL 

E 

CO 

CO 



CD 
LO 



LO 
CO- 
LO 





CO 


cz 


QQ 


"53 
E 


oo 

CD 

■o 


latio 


"O 
CZ 

CO 






CO 


time 


inci 


corr 


<C 

M 

o 



LO 




o o 

CD O 

=3 CD 

cr hz 

CD o 



LO 
1 — ' 
LO 




CO* 
LO 



LO 
CXI 
LO 



I 







cz 




o 






cz: 




IE? 


CD 






Z3 


CD 




czr 






CD 


CO 








cz 




CD 


.2 




CZ 




~z2 


CD 


~o 




ZD 


CD 




czr 






CD 


CO 




« trr CD 
Z3 

czr "3 

CD 



CD 
CO 

CO 

o 

CD 



«1 — CO 

cxj to 

^ CD 

E 
CO 



cr> 

CD- 
LO 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



PCT/EP99/03793 



6/31 




SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 



8/31 



PCT/EP99/03793 



CO 
CD 





=0.571 


=0.444 


=0.762 










+ 


+0.143 


-0.222 


-0.238 


CO 

~u 

j 


-0.143 


+0.333 


-0.429 


cnj 

+ 


-0.429 


-0.111 


+0.381 


+ 


CO 

co 

CNJ 






CD 

+ 


CD 

+ 


+ 


1 


CO 
CO 
CNJ 








CD 

I 


CD 

i 


CD 

i 


CNJ 

1 


+0.429 


+0.111 


-0.381 


CO 
1 


+0.143 


-0.333 


+0.429 


i 


-0.143 


+0.222 


+0.238 


< 


CNJ 
CO 

5 ii 

II a 


co 
S II 
II £ 
^ ST 


co ^> 

ss 'L, 

Co oo 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 PCT/EP99/03793 



9/31 



FIG. 10 



lm{l_C[0]} 
A f =1 




Re{L c [0]} 



Af=-1 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



10/31 



PCT/EP99/03793 



FIG. 1 1 




10log 10 (£ s /«o) [IB] 



o 


- -o 


CF 


|£ll<0.714 


X- 


X 


CF 


1^1 1 =0.8 


0- 


--0 


CF 


1^1 1 =0-9 


V- 


- -V 


CF 


1^1 1-1.0 






FD 


|^|<0.286 


X- 


X 


FD 


0.286 < | ^1 1<0.857 


o- 


—0 


FD 


0.857 < | ^ | < 1.429 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



11/31 



PCT/EP99/03793 



FIG. 12 




10log 10 (fsWo) [OB] 



o - 


-o 


CF 


|£l |<0.714 


X- — 


—X 


CF 


1^1 1 =0.8 


o- 


-o 


CF 


l$1 1 =0.9 


V- - 


-V 


CF 


1^1 1=1-0 






FD 


1<0.286 


X 


— X 


FD 


0.286 <|^ |<0.857 






FD 


0.857 < | <1. 429 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



12/31 



PCT/EP99/03793 




SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



13/31 



PCT/EP99/03793 



FIG. 14 



Pft 




10 12 14 16 
10log 10 (£ s //V 0 ) [dB] 



simulation, AWGN channel 



o 

X 

0 



- - theory, AWGN channel 

o CF: \^ |=0.0 0.7 

x CF: 1 =0.8 
o CF: | $1 1=0.9 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



14/31 



PCT/EP99/03793 



FIG. 15 




1 1 1 1 1 1 1 1 — I 

0 2 4 6 8 10 12 14 16 



10log 10 (£ s //V 0 ) [dB] 

simulation, channel profile 1 

theory, AWGN channel 

o o CF: 1 5i | =0.0,..., 0.6 
x xCF:|£i|=0.8 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 PCT/EP99/03793 



15/31 



FIG. 16 




simulation, AWGN channel 

theory, AWGN channel 

o o FD: | £i | =0.0, ...,0.2 

x x FD: | =0.3, .... 0.8 

o o FD: | $1 1 =0.9 1.4 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



16/31 



PCT/EP99/03793 



FIG. 17 




simulation, multipath channel 

theory, AWGN channel 

o o FD: | $1 1 =0.0. .... 0.2 

x x FD: | $1 1 =0.3, ...,0.8 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



17/31 



PCT/EP99/03793 



FIG. 18 




1 1 1 1 1 1 1 1 1 

0 2 4 6 8 10 12 14 16 



10log 10 (£ s //V 0 ) WB] 

simulation, (a prori) 

theory, (equiprobable) 

o o FD: | | =0.0 0.2 

x x FD: |^1 1 =0.3, ...,0.8 
o o FD: 1 =0.9 1.4 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 PCT/EP99/03793 

18/31 



FIG. 19 




10log 10 (£ s //Vo) [dB] 



CF: | =0.0 0.6 

CF: | =0.8 

o o D m =0 (ideal) 

v v Dm =1 

o o D m =2 

a a D m =3 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



19/31 



PCT/EP99/03793 



FIG. 20 




0 2 4 6 8 10 12 14 16 



10log 10 (£ s /A/ 0 ) [dB] 

FD: |4ll=0.0 0.2 

FD: 1 =0.4 0.8 

FD: |^|=1.3 

o o D m =0 (ideal) 

v v o m =1 

o o D m =2 

a a D m =3 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



20/31 



PCT/EP99/03793 



FIG. 21 




10log 10 (%Wo) [dB] 

FD: | ^ | =0.0 0.2 

FD: 1 =0.4 0.8 

FD: |^1 1 =1.3 

o o D m =0 (ideal) 

v v Dm =1 

o o D m =2 

a a D n =3 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



PCT/EP99/03793 



21/31 



FIG. 22 



Pft 




10" 4 -:e 



10-3-;EEHHEnb 



0 



2 



4 



6 



8 10 12 14 
10log 10 (£ s /A/ 0 ) [dB]- 



16 



o 

V 



o 

V 



FDPC: 1 =0.0, ...,0.8 
D m =0 (ideal) 
D m =2 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 PCT/EP99/03793 

22/31 



FIG. 23 




0 2 4 6 8 10 12 14 16 



10log 10 (£ s //V 0 ) [dB] 







FDPC: 1 =0.0, ...,0.8 


o 


O 


D m =0 (ideal) 


V 


V 


Dm=1 


0 


0 


Dm=2 


A 


A 


D m =3 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



PCT/EP99/03793 



23/31 




SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



24/31 



PCT/EP99/03793 




SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



25/31 



PCT/EP99/03793 



FIG. 27 




lower bound (theory) 

after iteration 1 

after iteration 2 

o o D m =0 (ideal) 

+ + D m =2 

x x D m =4 

* * Dpp=6 

o o D m =8 

v v D m =10 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



26/31 



PCT/EP99/03793 



FIG. 28 




10log 10 (£ s //Vo) »B] 



_ _ _ after iteration 1 

after iteration 2 

o o D m =0 (ideal) 

+ + D m =2 

x x D m =4 

* * Dpp=6 

o o D m =8 

v v D m =10 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



27/31 



PCT/EP99/03793 



FIG. 29 




8 



o 
+ 

X 

* 

O 
V 



o 
+ 

X 
* 

o 

V 



10 12 14 16 18 20 
10log 10 (£ s //Vo) [dB] 



after iteration 1 
after iteration 2 
D m =0 (ideal) 
D m =2 
D m =4 
D m =6 
D m =8 
D m =10 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



28/31 



PCT/EP99/03793 



FIG. 30 




-|0~6 H i i i i i i i i i 

0 2 4 6 8 10 12 14 16 18 20 

10log 10 (/y/Vo) [dB] 



after iteration 1 

after iteration 2 



o 


o 


D m =0 (ideal) 


+ 


+ 


Dm=2 


X 


X 


D m =4 


* 




D m =6 


0 


<> 


D m =8 


V 


V 


D m =10 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



29/31 



PCT/EP99/03793 



FIG. 31 




1 1 1 1 1 1 i 1 1 1 I I — H 

6 7 8 9 10 11 12 13 14 15 16 17 18 



10log 10 (£ s //V 0 ) [dB] 

ideal PER 

PER (coarse time estimate) 

PER (refined time estimate) 

o o 0 postfix symbols 

o o 1 postfix symbols 

* * 2 postfix symbols 

□ □ 3 postfix symbols 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



30/31 



PCT/EP99/03793 



FIG. 32 




6 7 8 9 10 11 12 13 14 15 16 17 18 



10log 10 (F s //V 0 ) [dB] 

ideal PER 

PER (coarse time estimate) 

PER (refined time estimate) 

o o coarse time: MMSE criterion 
* * coarse time: MC criterion 



SUBSTITUTE SHEET (RULE 26) 



WO 99/65180 



31/31 



PCT/EP99/03793 



FIG. 33 




i 1 1 1 1 1 i 1 1 1 1 1 — 

6 7 8 9 10 11 12 13 14 15 16 17 18 



10log 10 (£ s /A/ 0 ) [dB] 

ideal PER: time, frequency 

□■ — □ ideal PER: frequency (MMSE) 

PER (coarse time estimate) 

PER (refined time estimate) 

o o time direction 

* * frequency direction 



SUBSTITUTE SHEET (RULE 26) 



